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PREFACE 

The  following  thesis  is  a  presentation  of  the 
results  obtained  during  research  on  the  use  of  correlation 
techniques  for  automatic  machine  recognition  of  unformatted 

text.  This  thesis  should  be  read  as  the  seventh  volume  in  a 
series  of  eight  other  works  involving  such  research  by 
former  students  at  the  Air  Force  Institute  of  Technology 
(A FIT).  In  sequential  order  (along  with  their  author),  the 
other  theses  are: 

1.  PATTERN  RECOGNITION  BY  FOURIER  SERIES  TRANSFORMS 

C  .  H  .  R  a  d  o  y 

2.  THE  CLASSIFICATION  OF  VISUAL  IMAGES  BY  SPATIAL 
FILTERING 

O.  H.  Tallman 

3.  FOURIER  SPATIAL  FREOUENCY  COMPONENTS  AS  DESCRIPTORS 

FOR  RECOGNITION  OF  VARIABLE  FONT  TYPESCRIPT 

P.  W.  Kabler 

.  THE  DESIGN  OF  AN  OPTIMUM  ALPHANUMERIC  SYMBOL  SET  FOR 
COCKPIT  DISPLAYS 

L .  F .  Bush 

5.  SEGMENTATION  OF  TOUCHING  ENGLISH  LETTERS 

R.  E.  Bentkowski 

6.  MACHINE  SEGMENTATION  OF  UNFORMATTED  CHARACTERS 

R.  A.  Simmons 

8.  WHOLE  WORD  RECOGNITION  BASED  ON  LOW  FREQUENCY 
FOURIER  COMPLEX  AND  AMPLITUDE  SPECTRUMS 

M .  Oha i r 

These  other  theses  provide  additional  understanding 

as  to  why  certain  parameters  exist  and  why  correlation  was 

performed  with  Fourier  Transforms. 

During  the  course  of  this  research,  problems  were 
encountered  with  the  OCTEK  2000  and  the  Eclipse  system.  With 
the  Octek,  random  pixel  shifts  occurred  in  the  digitization 
of  pictures.  The  source  of  this  problem  was  still 
undetermined  at  the  ronclusion  of  this  thesis.  Problems  with 
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the  Eclipse,  however,  were  known  to  be 
Sharing  System)  operating  system.  Certai 
be  run  under  this  system. 

In  order  to  complete  this  thesi 
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ABSTRACT 

Th i s  thesis  presents  an  algorithm  tor  recognition  of 
letters  or  the  English  alphabet.  The  letters  may  be  single, 
isolated  letters;  or  concatenated,  groups  of  letters.  The 
algorithm  essentially  consists  of  a  two  pass  correlation  in 
which  the  size  and  energy  of  the  unknown  characters  are  set 
equal  to  a  known  template. 


INTRODUCTION 


S  i  n  .  ■  e  the  advent  >t  !i  !  g  h  -  s  i  «•  e  d  computers,  digital  image 
pro. -easing  has  become  a  topi  ot  great  interest  to  both 
a  .  a  d  e  m  1  a  and  i  .iiiust  rv  .  In  particular,  industry  has  t  o  >  ij  s  s  e  d 
much  attention  on  that  ot  machine  recognition  or  two-  and 
three-dimensional  images.  Although  much  hardware  and 
software  have  been  developed  during  the  last  tew  years  to 
support  such  e  t  torts,  modern  computers  still  cannot  locate 
and  identity  objects  within  a  cluttered  scene.  This  extends 
trom  locating  and  identifying  letters  on  a  page  to  t  o. 
tanks  in  the  desert 

This  thesis  addresses  the  former  problem.  Specifically, 
it  examines  the  problem  or  locating  and  identifying  printed 
Knglish  letters  on  a  page  ot  text. 


1 . 1  GENERAL  BACKGROUND 

Human  entry  of  data  into  digital  computers  is  often 
slow  and  erratic.  In  most  cases,  it  involves  the  conversion 
ot  printed  symbols  (useable  by  people)  into  electronic 
symbols  (useable  bv  computers).  A  secretary  usually  performs 
this  task  by  reading  the  input  data  and  then  toggling  a  set 
or  switches,  or  keys,  attached  to  a  keyboard  to  enter  the 
data  into  the  computer. 

Since  this  method  of  data  entry  is  totally  under  human 
control,  the  abilities  of  each  individual  operator  greatly 
affect  the  performance  and  operation  of  the  process.  The 


I 


physical  endurance  and  attention  span  o  t  e  a  e  h  person  limits 
the  number  o  t  hours  of  operation.  (Alter  several  hours  of 
work,  a  person  usually  becomes  bored  and  tired.)  In  addition 
to  human  fatigue,  the  process  is  also  affected  by  the  skills 
ol  each  individual  operator.  Thus,  speed  and  accuracy  varies 
accordingly  with  the  competence  of  each  person. 


r  n 

1  ight 

o  f 

the 

fact  tli  at  machines 

can  per 

form  at 

.iter 

speeds 

and 

tor 

longer  hours  than 

a  person 

,  there 

exists  a  need  for  machine  control  of  the  entire  data  entry 
process.  Automation  would  improve  efficiency  and  reliability 
and  reduce  the  total  amount  of  manpower  required  for  the 
task. 

!  .  1  A  BACKGROUND  on  data  entry 

At  present,  most  applications  require  that  data  (  i  . e  . 
hand-written  text  and  untormatted  text)  be  entered  manually 
by  human  operators  into  digital  computers  via  keyboards. 
This  requires  manv  hours  of  tedious  manual  labor.  An  ideal 
process  would  be  one  where  the  operator  merelv  presents  the 
computer  with  some  printed  material  and  it  automatically 
inputs  the  data.  For  ease,  the  prinLed  material  mav  be 
presented  visually  (pages  from  a  document)  or  ver  b  a  1  1  y 
('  s  p  o  k  e  n  nit  o  a  microphone). 

The  levelopment  of  such  alternatives  to  manual  keyboard 
entry  or  data  include  research  for  Optical  Character 
R  e  ,  o  g  n  i  t  l  o  n  (OCR)  systems  and  Voice  Recognition  systems.  For 
the  most  part,  it  appears  as  if  replacement  of  the  keyboard 


II.  THEORETICAL  DEVELOPMENT 


One  of  man's  greatest  desires  has  always  been  to  create 
a  machine  which  functions  in  human  like  manner.  Even  though 
he  still  does  not  fully  understand  the  mechanics  of  the 
ubiquitous  data  processing  operations  (which  seem  trivial 
and  commonplace  to  him)  that  he  routinely  performs,  he  will 
nevertheless  attempt  to  replicate  those  processes  in  a 
mechanical  device.  He  has  thusfar  produced  mechanical 
devices  which  can  record,  reproduce,  and  transmit  speech  and 
pictures.  Now,  he  must  develop  machines  which  can  read, 
understand  speech,  and  analyze  images. 

With  vision  as  the  subject,  several  researchers  have 
spent  much  time  and  effort  in  developing  models  for  the 
human  visual  system.  These  models  may  provide  the  design  of 
future  optical  recognition  machines. 

2 . 1  INITIAL  CONSIDERATIONS 

In  our  attempt  to  build  a  reading  machine,  we  are 
trying  to  imitate  the  human  process  ot  reading.  Therefore, 
the  machine  process  should  follow  that  of  the  human  process. 
However,  since  we  still  do  not  fully  understand  (or  even 
have  a  complete  theory  on)  how  the  human  visual  system 
processes  the  information  it  re  c  eives,  it  would  be  very 
difficult  to  design  a  recognition  device  which  employs  the 
s  a m e  variables  that  a  person  uses  to  describe  h  l  s 
e  n  v  1  r  u  n  m  e  n  t  . 


1.  DG  Eclipse  S/250 

(  1  6  -  b  1  t  minicomputer  with  256  kbytes  of  memmory  and 
-•  harddisk  drives) 

2.  DG  Nova  2 

(16-bit  minicomputer  with  64  kbytes  of  memmory) 

3.  Octek  2000  Image  Analyzer 

Dasher  60  5  3  terminal  t,  f  o  r  the  Eclipse) 

5.  Heathkit  terminal  (for  the  Nova) 

6.  Tektronix  *  6  3  2  video  hardcopy  unit 

7.  Cohu  vidicon  camera 
(model  6150-000) 

8.  Electrohome  display  monitor 


m. 


exists  or  can  be  easily  written  at  a  future  time. 

Such  programs  will  perform  the  following  tasks: 

a.  determine  the  angle  of  a  line  and  make  adjust¬ 

ments  for  it, 

b.  determine  the  height  of  a  character  and  adjust 

that  of  the  template. 

1 . 6  APPROACH 

Research  proceeded  from  the  idea  that  machine  reading 
should  be  based  upon  character  recognition,  i.e., 
identifying  each  letter  on  a  page.  It  was  also  determined 

that  recognition  should  be  performed  using  correlation 
techniques. 

In  order  to  simplify  the  recognition  task,  it  was 
decided  that  as  much  noise  as  possible  be  removed  from  the 

input  text.  Thus,  the  first  part  of  this  research  examined 

various  methods  for  the  automatic  reduction  and  elimination 
of  noise. 

1 . 7  MATERIAL  AND  EQUIPMENT 

The  Signal  Processing  Laboratory  at  the  Air  Force 
Institute  of  Technology  (AFIT)  provided  all  the  necessary 
material  and  equipment  for  this  research.  Guidance  was  given 

bv  AFIT  staff  and  faculty. 


Most 

of  the  work 

was  performed 

o  n 

a  Data  General 

(  DC  ) 

computer 

system  which 

consisted  of 

the 

following  : 

1 
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uncorrelated,  that  is,  each  digitization  of  the  same 
text  will  produce  a  random  noise  pattern.  If  the 
noise  pattern  is  constant,  then  a  simple  filter  can 

mask  it  out. 

3.  Random  noise  will  not  alter  a  character's  shape, 
i.e.,  stray  marks  will  not  change  one  letter  into 
another.  Usually,  stray  marks  are  much  lighter  in 
color  than  that  of  the  letter  and  can  be  detected  as 

noise.  For  example,  stray  pencil  marks  will  not 
change  a  " c "  into  an  "o".  However,  if  such  marks 

are  darker  than  the  letter,  it  will  be  recognized  as 

an  " o " . 

*.  The  input  text  is  properly  lighted,  such  that, 
•  digitization  of  the  text  will  produce  a  bi-modal 

distribution  of  greylevels.  Much  of  the  background 

will  then  be  centered  around  white,  and  much  of  the 

print  will  be  centered  around  black. 

5.  It  is  assumed  that  a  mechanical  device  exists 

which  can  automatically  control  the  focussing  of  the 

digitizing  camera.  Such  devices  are  available 
c  omme  rcially. 


The 

lines  of  printed 

text 

are  placed 

such 

that 

each 

line  is  horizontal 

and 

not  at  an  an 

g  1  e  . 

Thus 

i  t  i 

s  assumed  that  the 

text 

is  properly 

placed 

f  o 

digitization  and  not  upside  down. 

assumption  is 

1  3 


The  final 


that 


certain 


sof  tware 


which  implemented  most  of  the  functions. 

Due  to  time  constraints,  this  thesis  could  not  examine 
every  aspect  of  the  recognition  problem  associated  with 
every  character  shape.  This  limited  the  study  only  to  typed 
letters  of  the  English  alphabet..  This  excluded  stylized 
print  (Script,  Gothic,  etc.),  numerics  (1,  23,  58,  etc.), 

and  special  characters  (*,  &  ,  %  ,  it  ,  etc.).  In  particular, 

this  thesis  did  not  deal  with  hand  written  text  or  print 
which  has  been  rotated.  For  most  of  the  testing,  a  book  was 
randomly  taken  out  of  the  library  and  pages  from  it  were 
digitized. 

The  only  constraint  set  upon  the  shape  of  the  letters 
was  that  each  character  must  be  recognizable  by  a  person 
without  reference  to  ..he  word  itself.  Thus,  recognition  of  a 
letter  must  be  based  solely  on  the  shape  of  the  letter  and 
not  on  what  it  should  be  in  a  specified  word. 

1 . 5  assumptions 

Seven  assumptions  were  made  at  the  start  of  this 
research  in  order  to  simplify  the  problem: 

1.  The  average  width  of  a  character  is  a  given  para¬ 
meter.  For  the  moment,  this  is  a  valid  assumption 


bee 

a  u  s  e 

i  t 

a  1  lows 

us  to 

search  for  a 

s  o  1  u  t 

ion. 

0  n 

any 

page 

o  f 

text. 

there 

is  a  high  probabili 

ty 

that 

one 

can 

find 

at  least  one 

one-lettered 

word. 

I  f 

one 

i  s 

found 

,  then  its 

width 

will  be  used. 

2.  The  noise  produced  by  the  digitizer  is 


classified  as  spatial  noise.  Noise  incurred  during  the 
digitization  process  include  video  snow  and  nonlinear  image 
distortions.  Such  noise  are  classified  as  temporal  noise. 
Noise  from  either  process  tends  to  alter  enough  of  the  shape 
of  the  letters  to  cause  errors  in  recognition. 

Other  problems  for  reading  machines  come  from  varying 
font  styles.  All  text  do  not  come  from  the  same  printing 
press.  Even  within  the  same  book  one  can  find  multiple  font 
styles.  (For  instance,  the  title  of  a  book  and  its  text  are 
usually  printed  in  different  styles). 

In  addition  to  varying  the  character  shape,  modern 
electronic  printing  with  its  proportional  spacing  causes  the 
virtual  concatenation  of  letters  within  words.  This  now 
creates  the  problem  of  finding  or  separating  the  letter  to 
be  recognized  ("segmentation"  problem). 

This  thesis  examines  those  problems  associated  with  the 
proccess  of  recognizing  letters.  Specifically,  it  examines 

1.  noise  from  the  digitization  and  printing  process, 
and 

2.  character  recognition  based  upon  correlation 
techniques  . 

1 . A  SCOPE 

The  goal  of  this  research  was  to  start  development  of 
an  algorithm  which  would  give  computers  the  ability  to 
recognize  the  letters  from  any  page  of  text  without  the 
intervention  of  a  human  operator.  The  primary  output 
consisted  of  the  algorithm  and  a  series  of  computer  programs 


1.3  STATEMENT  OF  PROBLEM 


To  date,  reading  machines  based  upon  Curve  Tracing  and 
Template  Matching  techniques  do  not  work  well  unless  the 
input  text  is  very  precisely  controlled  in  shape  and 
location.  The  problem  lies  in  the  fact  that  such  techniques 
fail  when  noise  is  encountered  in  the  print  which  causes  the 
input  font  style  (type  style)  to  vary  from  that  of  the 
standard  set.  In  addition,  current  techniques  do  not  allow 
for  the  identification  ofletters  which  do  not  come  from  the 
standard  set. 

The  performance  of  current  reading  machines  is  optimum 
only  when  the  print  on  the  input  text  is  "perfect". 
"Perfect"  means  the  following  conditions: 


1  . 

the 

ink  is  of  a 

uni  form 

black  color, 

2  . 

the 

paper  is  of 

a  uniform  white  color. 

3  . 

the 

letters  are 

all  from 

the  same  type  style 

4  . 

the 

letters  do 

not  touch 

each  other,  and 

5  . 

the 

letters  are 

all  well 

f  ormed  . 

Under  these  conditions,  a  reading  machine  will  not  have  any 
difficulty  in  recognizing  the  letters  on  a  page  of  printed 
text. 

Most  text,  however,  does  not  come  in  perfect  condition. 
One  source  of  problems  is  that  of  noise  incurred  during  the 
printing  process  and  the  digitization  process.  Noise 
incurred  during  the  printing  process  include  splotches  of 
ink,  strav  marks,  and  breaks  in  the  letters.  Such  noise  are 


1  .  2  D  TECHNOLOGICAL  LIMITATIONS 

Using  today's  more  advanced  technology,  computer  vision 
can  still  at  best  be  compared  to  that  of  a  blind  man  probing 
around  his  environment  with  a  cane.  With  respect  to  reading 
unformatted  text,  the  scene  becomes  that  of  the  blind  man 
probing  around  a  page  of  raised  letters  with  a  pin.  Both 
the  man  and  the  machine  cannot  examine  an  entire  page  at 
once  to  locate  the  print.  They  must  search  around  the  page 
until  a  line  of  print  is  found.  Then,  they  must  find  the 
beginning  of  the  line  and  determine  the  first  letter  of  the 
first  word. 

Restricted  to  using  only  a  pin,  finding  and  then 
recognizing  letters  becomes  a  very  difficult  task.  for  the 
man.  (Allowing  him  to  use  his  finger,  or  hand,  would  only 
mean  programming  the  computer  with  more  pointers  and 
performing  more  computations.)  Although  it  is  a  difficult 
task,  it  is  not  an  impossible  task  for  the  man.  Thus, 
building  a  machine  that  reads  should  also  not  be  an 
impossible  task. 

Today's  technology  can  provide  all  the  computational 
tools  required  to  perform  the  necessary  calculations.  In 
conjunction  with  existing  mass  storage  systems  and  high 
resolution  digitizers,  all  the  hardware  needed  to  build  a 
reading  machine  exists.  However,  complete  software  required 


to  run  the  hardware  does  not  exist 


The  Kurweil  Computer  Products  company,  producers  of  the 
Kurweil  Data  Entry  Machine,  keeps  the  designs  for  their 
machine  a  closely  guarded  secret.  The  techniques  it  employs 
are  not  made  known  to  the  public. 

Although  they  are  sold  commercially,  most  machines  do 
not  perform  as  well  as  expected.  Highly  varying  font  styles 
and  nonuniform  printing  cause  intolerable  error  rates.  A 
human  operator  is  then  required  to  correct  the  mistakes. 
Thus,  instead  of  c  t  rating  at  machine  speeds,  such  machines 
operate  at  rates  slightly  higher  than  that  of  a  secretary. 

1 . 2C  AFIT  OCR 

At  the  Air  Force  Institute  of  Technology  (AFIT),  one 
technique  for  OCR  employs  a  Fourier  Transform  correlation 
scheme  for  its  recognition  process.  Developed  by  R.  Simmons 
in  1981,  the  algorithm  reconstructed  the  the  magnitude  of 
the  Fourier  Transform  of  a  template  character  with  the  phase 
of  the  Fourier  Transform  of  a  string  of  characters  (which 
contains  the  template  character).  Shapes  which  did  not  have 
the  magnitude  of  the  template  were  attenuated  in  the  visual 
domain.  Thus,  a  peak  detector  could  be  used  to  locate  the 
most  probable  position  of  a  character. 

This  algorithm  seemed  to  work  quite  well  with  "large" 
letters,  but  failed  when  presented  with  letters  that  were 
"smaller"  (in  comparison  with  the  Fourier  window  size). 


standard  set  for  recognition. 


1  .  2  A  APPLICATIONS  FOR  READING  MACHINES 


A  s 

stated  in  the  previous  sections 

,  the 

p  r  i  ma  r  y 

use  of 

a  r  e  a  d  l 

n  8 

machine  is  for 

au  t  oma  tic 

entry 

of  data 

into 

computers  . 

This  includes 

use  for 

text 

storage 
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informat 

i  o  n 

processing.  One 

specific  area  of 

appl icat 

ion  in 

in  format 

ion 

processing 

involves 

automatic  computer 

translation  of  foreign  text.  As  a  front-end  processor,  a 
reading  machine  can  automatically  identify  and  input  foreign 
text  into  a  computer  which  can  then  translate  it. 

A  very  well  known  application  area  for  reading  machines 
involves  the  Post  Office  and  sorting  mail.  Such  a  device  can 
identify  the  address  on  a  letter  and  classify  each  with 
respect  to  their  zip  code. 

In  addition  to  easing  the  data  entry  task,  reading 
machines  can  also  be  used  as  aids  to  those  visually 
impaired.  Connected  to  a  voice  synthesizer,  such  machines 
can  literally  read  to  persons  who  cannot  see. 

1  .  2  B  OPERATIONAL  READING  MACHINES 

At  the  present,  there  are  several  reading  machines 
available  commercially.  They  include  the  RCA  Multi-Font 
Reading  Machine,  the  IBM  Reading  Machine,,  and  the  Kurzweil 
Data  Entry  Machine. 

The  RCA  machine  uses  an  optical  correlation  technique 
in  its  recognition  scheme.  A  "picture"  of  an  unknown 
character  is  projected  onto  a  transparency  containing 
negative  and  positve  replicas  of  all  the  possible  characters 


3.  classify  Che  objects  within  the  visual  image 
(e.g.  determine  if  the  object  is  a  tank). 

Specifically  for  reading  machines,  a  video  digitizer 
converts  the  pages  of  "  input  text  into  digital  images.  A 
computer  or  microprocessor  then  analyzes  each  digital  image 
for  characters. 

It  is  important  to  note  that  the  digital  image,  or 
digitized  picture,  is  a  bit  map  of  the  light  intensities 
reflected  from  the  page.  Therefore,  the  electronic  image  is 
only  a  "photograph"  of  the  text.  It  contains  only  greylevel 
information  and  does  not  contain  ASCII  values  for  the 
letters  on  the  page.  A  computer  or  microprocessor  is 
required  to  process  the  greylevel  image  to  find  and  identify 
the  letters.  Usually,  features  extracted  from  the  electronic 
images  are  compared  against  some  set  of  templates  and 
thresholded  to  identify  the  letter. 

1.1C  background  on  processing  techniques 

Traditional  methods  for  processing  characters  have 
included  Curve  Tracing  and  Template  Matching  techniques. 
With  each,  a  prototype  (or  token)  is  generated  with  some 
specific  attributes  of  a  given  letter.  Prototypes  for  each 
of  the  given  letters  then  form  a  template  set  (or  standard 
set).  Accordingly,  each  input  image  is  compared  against 
those  in  the  template  set  and  categorized  as  being  similar 
to  one  or  more  prototypes;  else  sufficiently  dissimilar  so 
as  to  be  rejected  as  an  allowed  character. 


4 


will  come  from  OCR  technology  in  the  form  of  a  reading 
machine. 

One  major  disadvantage  of  a  v o i c e - o p e r a t e d  input  device 
is  that  it,  too,  operates  under  human  control.  (The  operator 
has  only  exchanged  the  use  of  one  set  of  muscles  for  those 
of  another.)  Thus,  like  keyboard  entry,  such  a  system  would 
be  susceptible  to  the  physical  limitations  of  the  individual 
operator.  However,  for  reading  machines,  control  is  under 
machine  guidance  and  operation  may  proceed  at  machine  speeds 
with  limitations  due  only  to  mechanical  performance. 

1.1B  BACKGROUND  ON  reading  machines 

In  this  thesis,  a  reading  machine  is  defined  as  any 
optical/electrical/mechanical  device  designed  to  recognize 
characters  from  a  given  set  of  letters  for  the  purpose  of 
automatic  input  of  data  into  digital  computers.  (See  setion 
1.2A  for  other  applications.)  Currently,  most  reading 
machines  are  optical  systems  based  upon  Optical  Character 
Recognition  technology. 

Since  reading  machines  deal  with  identifying  letters, 
they  are  classified  as  image  processing  devices.  The  basic 
operation  of  all  such  devices  consists  of  the  following 
steps  : 

1.  convert  the  visual  image  into  a  digital  image 
(e.g.  a  greyscale  version  of  the  visual  image) 

2.  store  and  pre-process  the  electronic  image 
(e.g.  edge-enhance  the  image) 


I 


We  do  not  know  how  people  read;  but,  we  do  know  that 
people  can  read.  The  human  reading  process  involves  input 
through  the  visual  system  and  processing  in  the  brain.  We 
must  therefore  consider  the  operations  of  these  systems  and 
use  them  as  guides  in  building  a  reading  machine. 

2 . 1 A  GUIDELINES  FOR  recognition 

If  the  human  recognition  system  exemplifies  the  best 
recognition  system,  then  our  mechanized  model  should  operate 
on  a  parallel  level. 

Under  certain  circumstances,  our  human  recognition 
svstem  tails  and  produces  erroneous  outputs.  Under  other 
conditions,  our  human  system  produces  the  proper  results. 
Thus,  given  the  same  conditions,  a  mechanical  device  should 
see  and  interpret  things  only  as  we  see  them. 

Remembering  that  the  reading  machine  is  a  mechanical 
device,  we  should  expect  it  to  outperform  the  biological 
machine  in  certain  areas.  This  comes  about  from  the  fact 
that  many  of  the  parameters  in  a  mechanical  device  are 
easily  controlled.  Hence,  problems  which  may  appear 
difficult,  or  even  impossible  to  the  human  system  should 
present  no  problem  to  the  mechanical  one.  (For  example,  a 
person  may  distinguish  between  only  approximately  100 
greylevels,  but  a  machine  may  be  programmed  to  distinguish 
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more.  ) 


2. IB  STANDARDS  FOR  OPERATION 


It  a  human  operator  can  isolate  and  identify  any  letter 
on  a  page  of  text,  then  a  reading  machine  must  also  be  able 
to  perform  the  same  task.  In  addition,  the  machine  must 
concur  with  the  letter  that  the  human  operator  recognized 
and  do  so  consistently. 

Since  the  mechanical  device  is  to  replace  a  human 
operator,  then  the  machine  must  perform  more  efficiently 
than  a  person.  It  must  have  a  higher  rate  of  operation  and  a 
lower  error  rate. 

2.2  READING 

We  do  not  know  exactly  what  sequence  of  mental 
activities  occurs  within  a  person  during  the  reading  process 
that  gives  him  the  ability  to  recognize  letters  of  varying 
sizes  and  styles.  (A  mere  child  of  five  may  be  taught  to 
recognize  the  26  letters  of  the  English  alphabet,  yet,  a 
computer  cannot  perform  the  same  task.)  In  mechanizing  the 
reading  process,  we  must  consider  the  method  of  recognition 
and  its  implementation. 

2  .  2  A  WORDS  OR  LETTERS 

Research  in  developing  reading  machines  (at  AFIT)  has 
recently  taken  on  two  approaches.  One  view  holds  the 
position  of  building  reading  machines  based  upon  word 
recognition,  i . e . ,  locate  the  words  on  a  page  and  try  to 
recognize  them.  The  other  and  more  traditional  approach 


bases  recognition  upon  the  letters  in  a  word. 

The  first  approach  is  much  simpler  in  that  one  does  not 
have  to  deal  with  finding  characters  in  a  concatenated 
sequence  of  letters.  In  most  text,  the  words  are  fairly 
separated  from  each  other;  while  the  letters  are  spaced 
fairly  close  to  each  other.  The  biological  reading  process 
seems  to  tolerate  closely  spaced  words  as  long  as  they  do 
not  physically  touch  each  other.  Letters,  on  the  other  hand, 
may  touch,  but  not  overlap  more  than  about  10'/.  of  their 
width  (2). 

At  the  present,  we  cannot  tell  which  approach  will  lead 
to  the  building  of  the  first  successful  reading  machine.  In 
this  thesis,  we  will  take  the  latter  approach.  Advanced 
reading  skills  may  perhaps  take  place  at  the  word  level, 
but,  a  human  reader  can  always  recognize  the  individual 
letters  on  a  page  of  text.  (For  this  reason  alone,  research 
was  devoted  to  character  recognition.) 


2  .  2 B  BUILDING  A  TEMPLATE  SET 

N'o  two  objects  in  this  physical  universe  are  ever 

exactly  alike.  Thus,  the  human  recognition  system  acts 
accordingly  to  generate  a  set  of  standards  or  references  for 
interpreting  and  classifying  those  inputs  from  its  visual 
e  ironment. 

From  the  fact  that,  at  any  instant  of  time,  the  human 

system  can  intelligently  ignore  all  extraneous  details  which 


are  unimportant  to  the  recognition  process,  some  researchers 


have  suggested  that  the  biological  recognition  system 
generates  a  rather  loose  set  of  standards  which  is  ill- 
defined  and  subject  to  refinement  and  adjustment.  The  exact 
process  whereby  it  generates  and  modifies  those  standards  is 
still  unknown  to  man. 

From  a  totally  different  point  of  view,  some  believe 
that  every  visual  experience  is  stored  somehow  in  toto. 
Thus,  the  human  recognition  system  compares  images  from  an 
external  scene  with  those  of  previous  impressions  or 
experiences.  There  are  still  others  who  believe  that  the 
answer  lies  somewhere  in-between  those  two  extremes.  In  any 
case,  we  still  do  aot  have  the  total  answer. 

2  .  2  C  STORAGE  BY  COURIER  TRANSFORMS 

Based  upon  biological  and  physiological  information 
about  the  human  visual  system,  Kabrisky  (1)  has  developed  a 
theoretical  model  for  the  human  recognition  system.  In  this 
model,  all  "remembered"  patterns  (experiences)  are  stored 
away  topologically  (unaltered  under  a  homeomo r ph i sm )  .  The 
system  then  performs  its  recognition  by  correlating,  on  a 
po i n t - b y - po i n t  basis,  a  given  input  pattern  with  all  the 
stored  patterns  and  selecting  the  pattern  which  correlates 
the  best. 

As  pointed  out  by  Kabrisky,  the  amount  of  memory 
required  by  this  model  to  store  all  of  the  patterns  which  an 
average  person  can  remember  would  seem  to  be  prohibitively 
large.  One  solution  to  this  memory  problem  would  involve 


addition,  the  size  of  print  can  be  quickly  enlarged  or 
reduced.  With  the  proper  manipulations,  input  data  can  be 
transformed  into  a  form  more  amenable  to  the  recognition 
s  c  heme  . 

2  .  3 A  NOISE  REDUCTION 

Before  we  begin  to  read  a  book,  we  usually  turn  on  a 
light.  This  allows  us  to  see  the  printed  text.  In  much  the 
same  way,  we  must  first  ensure  that  the  reading  machine  can 
"see"  the  text.  Proper  illumination  and  placement  of  the 
text  for  digitization  provides  one  means  of  this. 

Once  digitized,  the  greyscale  distribution  can  be 
modified  for  maximum  contrast  between  print  and  background. 
Assuming  a  bimodal  distribution  of  color,  a  greyscale  image 
can  be  thresholded  to  form  a  binary  (black  on  white)  image. 
Binary  images  are  more  tractable  for  computer  processing. 

2 .  3  B  LOCATING  THE  LETTERS 

After  an  electronic  image  has  been  formed,  the  computer 
mu9t  read  the  print  on  the  page.  Following  the  human 
process,  the  computer  must  start  at  the  beginning  of  the 
text;  reading  the  first  word  of  the  first  line.  Here, 
instead  of  recognizing  the  entire  word,  the  process  is 
brought  to  a  lower  level  of  recognizing  the  letters  in  the 
first  word. 

In  order  to  read  the  letter,  the  machine  must  first 
find  the  location  of  the  words.  This  requires  finding  the 
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2.3C  RECOGNITION  OF  LETTERS 


Based  upon  Kabrisky's  model,  recognition  should  come 
from  correlating  the  Fourier  transform  of  an  input  charac¬ 
ter  with  those  of  a  template  set.  Also,  such  correlations 
should  take  place  in  the  Fourier  domain  with  spatial 
f i ltering. 

If  the  human  system  uses  an  loosely  defined  template 
set  for  recognition,  then  the  mechanical  model  should  be 
able  to  modify  its  templates.  Futhermore,  since  correlation 
is  similar  to  template  matching,  it  would  be  desirable  to 
have  both  the  template  shape  and  input  shape  as  similar  as 
possible.  This  implies  adjustment  of  the  height  and  width  of 
each  template  letter.  The  template  height  can  be  set  to  the 
input  height  by  measuring  that  of  the  input  shape. 
Adjustment  of  character  width  means  adjusting  the  thickness 
of  each  character. 

Adjusting  the  character  width  will  help  correlation  in 
that  this  will  affect  the  low  order  harmonics  of  a  Fourier 
Transform.  (A  thick  character  will  have  more  energy  content 
in  the  low  order  harmonics  than  a  thinner  version  of  the 
same  letter.)  Either  the  template  letters  can  be  thickened 
to  match  that  of  the  input,  or,  the  input  can  be  thinned  to 
match  that  of  the  template.  The  reasons  for  choosing  a 
thinning  algorithm  is  two-fold: 

1.  Thinning  offers  a  more  controllable  operation.  It  is 

easier  to  thin  a  letter  to  a  known  pixel  size  than  to 


thicken  a  letter  to  a  varying  pixel  size.  (Some  letters 
are  not  of  the  same  uniform  thickness  in  certain 
fonts.  ) 

2.  It  is  safer  to  modify  the  input  character  thickness 
than  vice  versa.  With  a  template  thicker  than  any  input 
letter,  the  possibility  exists  of  correlation  with 
nearly  everything  in  an  input  scene.  (See  Chapter  III.) 

With  the  height  and  width  adjusted  as  much  as  possible, 
the  remaining  parameter  to  adjust  is  the  energy.  If  an  input 
character  and  a  template  contain  the  same  energy,  then 
correlation  will  produce  a  maximum  peak  only  if  the  two 
figures  truely  have  the  same  shape.  Other  shapes  which  may 
be  similar  will  produce  lower  peaks. 


III.  SOFTWARE  DEVELOPMENT 

The  following  chapter  formally  presents  the  recognition 
algorithm  examined  in  this  thesis  research.  Although  a  total 
algorithm  was  described  in  Chapter  2,  only  a  portion  of 
those  functions  were  fully  implemented  in  software.  Rather 
than  devote  much  of  the  time  to  coding,  it  was  decided  that 
time  was  better  spent  on  testing  the  validity  of  the 
recognition  algorithm. 

This  chapter  is  divided  into  two  sections.  The  first 
section  introduces  the  overall  processing  procedure  of  the 
character  recognition  algorithm  implemented  in  this 
thesis.  The  second  section  briefly  expands  upon  the  steps  of 
the  processing  procedure  and  outlines  the  actual  operational 
process  of  the  recognition  scheme.  (Chapters  A  and  5 
describe  in  more  detail  the  operation  of  the  software 
implementation. ) 

3 . 1  OVERALL  PROCESSING  PROCEDURE 

The  entire  processing  procedure  of  the  recognition 
algorithm  can  essentially  be  broken  down  into  three  stages: 
the  digitization  stage,  the  pre-processing  stage,  and  the 
recognition  stage.  Each  stage  performs  its  operations 
independently  of  the  others  with  the  output  of  one  stage 
feeding  directly  into  that  of  the  next  one  (see  Figure  3.1). 
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FIG.  3.1  OVERALL  PROCESSING  PROCEDURE 


The  Digitization  stage  performs  the  actual  conversion 
of  analog  video  text  into  digitized  greyscale  text.  Most  of 
this  stage  consisted  of  hardware  which  had  already  been 
established  by  former  thesis  students  in  the  Digital  Signal 
Processing  Lab  at  AFIT. 

The  Pre-Processing  stage  takes  digitized  greyscale  text 
and  thresholds  the  greylevels  to  produce  binary  text.  In 
addition  to  this,  it  determines  the  locations  of  the  lines 
of  print  on  the  text  and  the  position  of  the  words  within 
the  lines. 

The  Recognition  stage  recursively  correlates  input 
characters  with  template  characters  in  order  to  identify  the 
letters.  At  each  iteration,  the  characters  are  energy 
normalized  to  produce  proper  correlation  peaks. 

The  reason  for  this  specific  division  of  stages  was 
twofold.  As  already  mentioned,  the  first  stage  was  already 
implemented.  There  was  no  need  to  duplicate  the  work  of 
others.  Also,  with  this  arrangement,  time  could  be 
appropriately  proportioned  to  each  of  the  other  stages.  If 
the  Pre-Processing  stage  was  not  fully  implemented  by  the 
alloted  time,  then  its  completion  would  not  affect  that  of 
the  Recognition  stage. 

3.2  THE  RECOGNITION  PROCESS 

The  recognition  algorithm  examined  in  this  thesis 
consisted  basically  of  correlating  an  unknown  word  with  a 
template  character.  (It  should  be  pointed  out  here  that  in 
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this  section,  and  those  that  follow,  the  term  "word"  refers 
to  the  output  of  the  Pre-Processing  stage.  A  word  can  be 
either  an  individual  character  or  a  string  of  concatenated 
letters.  For  an  explanation  of  this,  see  section  ‘4.2.)  The 
energy  of  the  template  which  produced  a  high  correlation 
peak  and  the  energy  in  the  area  of  the  word  where  the  peak 
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second  correlation  pass  would  indicate  the  most  probable 
identity  of  the  word.  (Maximum  in  the  sense  that  if  two 
objects  have  the  same  shape,  then  normalizing  their  energies 
would  produce  a  correlation  peak  of  one.) 

More  explicitly,  the  correlation  process  consisted  of 
the  following  steps: 

1.  Individually  correlate  the  characters  of  a  template 
set  with  an  unknown  word. 

2.  Determine  the  location  where  each  template 

correlated  well. 

3.  Center  a  window  around  the  correlation  peak  area 
and  set  the  energy  equal  to  that  of  a  template. 

■4.  Re-correlate  the  energy  modified  word  with  the 
template  . 

5.  Repeat  the  energy  modification  and  re-correlation 
process  for  each  of  the  templates  determined  in  the 
first  pass. 

6.  Compare  the  correlation  peak  values  of  the  second 
pass  and  choose  the  highest  as  the  best  answer. 
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IV.  THE  DIGITIZATION  STAGE 


This  chapter  discusses  the  digitization  stage.  As 

mentioned  before,  this  level  of  the  algorithm  had  already 

been  implemented  in  hardware  by  other  students.  Many  o  t  the 
programs  for  running  the  digitization  hardware  were  readily 
available.  Instead  of  describing  the  details  of  those 
written  programs,  this  chapter  highlights  those  features  of 
the  digitization  stage  which  are  important  to  the 
understanding  of  the  software  for  the  next  two  stages. 

..1  ESTABLISHED  WORK 

Most  of  the  original  software  :or  this  stage  was 

written  b  v  Robin  E.  Simmons.  As  a  thesis  student,  he  set  the 
toundation  and  established  the  useful  software  systems  for 
future  work  in  the  Digital  Signal  Processing  Lab.  The 

programs  he  wrote  included: 

1.  L'NPACK.fR 

2.  REPACK.FR 

3 .  VTOC . FR 
.  I  T  0  C  .  F  R 

These  were  svstem  dependent  Fortran  subroutines  which 
allowed  for  the  manipulation  of  picture  elements  (pixels) 

within  a  digitized  picture.  (See  section  below  for 

explanation. ) 

Since  the  software  alreadv  existed,  this  stage  was  used 
as  a  black  box  which  converted  analog  video  pictures  into 
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The  fourth  step  (FWORD.FR)  extra  c  ts  one  word  from  the 
line  ot  printed  text.  The  word  chosen  is  up  to  the  • .  r  e l  .  It 
copies  this  word  onto  a  blank  background  for  the  next  step. 

The  tilth  and  final  step  (FCHAR.FR)  further  divides  the 
word  it  it  is  possible.  The  third  step  only  examines 
vertical  blank  lines  between  words.  If  an  angular  line 
exists,  then  it  would  be  missed.  The  fifth  stage  accounts 
tor  that.  It  performs  an  operation  much  like  that  of  finding 
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Figure  5.7  EXTRACTION  OF  WORDS 


5.3  LINE/WORD  LOCATING  ALGORITHM 


The  input  to  the  LINE/WORD  locating  process  is  the 
noise  reduced  text.  From  this,  the  algorithm  determines  the 
locations  of  the  lines  of  print  and  proceeds  to  extract  the 
words  from  the  lines.  (In  this  section  and  those  which 
follow,  the  term  "word"  refers  to  the  output  of  the 
LINE/WORD  locating  process.  Thus,  a  word  may  be  an 
individual  character  or  a  string  of  concatenated  letters.) 

This  stage  performs  its  operations  in  five  steps  (see 
Figure  5.7).  The  first  step  (LINES. FR)  determines  the 
locations  of  the  lines  of  print.  It  does  so  by  calculating  a 
horizontal  histogram  of  the  text.  Assuming  that  the  text  is 
properly  placed,  there  should  exist  a  series  of  blank 


v  e  r  t  i 

cal 

lines  between  each 

line  of 

print  . 

.  These  b 1 ank 

lines 

would 

demarcate  the  location 

of  the 

lines 

of  print. 

The 

second  step  (FLINE 

.FR)  extracts 

one  line  of 

print 

from 

the 

text.  The  choice 

of  which 

line 

is  up  to  the 

user. 

It  copies 

the  specified  line 

onto  a 

blank 

background  t" 

or  the 

next 

step 

The  third  step  (WORDS . FR )  calculates  a  vertical 
histogram  of  the  chosen  line.  If  a  straight  vertical  line 
can  be  drawn  between  each  character,  then  that  character  is 
counted  as  one  word.  As  mentioned  above,  if  the  letters  are 
concatenated,  then  the  string  of  letters  is  considered  to  be 
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one  word. 


CLEAN  takes  a  16  level  greyscale  image  and  threshold- 
averages  it  to  a  binary  image.  It  performs  this  operation 
under  the  assumption  that  the  greylevels  have  a  bimodal 
distribution.  This  routine  first  calculates  a  histogram  of 
the  greylevels.  It  then  determines  the  two  peaks 
corresponding  to  character  color  and  background  color  (see 
Figure  5  .  .  From  the  peaks,  it  determines  the  first  valley 
after  the  first  peak.  It  uses  this  valley  point  to  begin 
thresholding  of  the  greyscale  image.  It  thresholds  up  to  the 
tenth  level  (the  start  of  white  on  the  scale).  After  each 
threshold  level,  the  routine  averages  the  clipped  values  to 
produce  a  t h r e s h o 1 d - a v e r a g e d  value.  This  final  value  is  then 
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The  noise  reduction  algorithm  developed  to  binarize  the 
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process  on  the  binarized  text-  It  examines  the  page  of  text, 

determines  the  locations  of  the  lines  of  print,  and  extracts  J 

a  word  from  one  of  the  lines. 

5.2  NOISE  REDUCTION  ALGORITHM  j 

In  order  to  understand  why  a  noise  reduction  routine  is 
necessary,  let  us  first  examine  the  output  of  the 

Digitization  stage.  Figure  5.1  shows  the  original  analog  j 

text.  Figure  5.2  shows  the  raw  digitized  greyscale  text. 

(The  text  had  been  properly  positioned  for  digitization  and 

illuminated  with  a  120  watt  lamp.)  j 

4 

The  greyscale  version  of  the  print  contains  unwanted 
temporal  noise  from  the  digitizer.  This  shows  up  as  the 

sprinkling  of  black  and  white  dots  on  the  text.  In  addition  | 

to  generating  this  type  of  noise,  the  quantization  error  of 

the  device  causes  the  distortion  of  greylevels  within  each 

character.  The  greylevel  value  of  each  pixel  within  a  letter 

can  vary  from  one  to  two  levels  of  black.  Clearly,  the 

contrast  between  the  characters  and  the  background  is 

greater  in  the  original  picture. 

Conversion  of  the  greyscale  text  into  binary  text  would 
allow  for  the  use  of  readily  available  digital  image 
processing  techniques.  (These  include  line  thinning  and 
curve  tracing  algorithms.)  Greyscale  text  would  limit 
processing  to  techniques  which  do  not  operate  on  the 


individual  pixels  of  a  picture. 


THE  PRE-PROCESSING  STAGE 


V  . 


The  following  chapter  discusses  the  software 
implementation  of  the  second  level  of  the  recognition 
algorithm.  (See  Appendix  A  for  complete  source  code 
listings.)  At  this  point,  analog  text  has  already  been 
digitized  and  is  contained  in  a  video  file.  The  Pre¬ 
processing  stage  will  strip  off  words  from  the  text  and 
present  them  to  the  next  stage  for  identification. 

The  first  section  of  this  chapter  describes  the  overall 
processing  function  performed  by  the  Pre-Processing  stage. 
The  second  and  third  sections  expand  upon  this  overview  and 
provide  more  details  as  to  the  actual  operation  of  each  of 
the  individual  programs  written  for  implementation. 

5 . 1  OVERALL  PROCESS 

The  Pre-Processing  stage  accepts  as  input  the  digitized 
greyscale  text  from  the  Digitization  stage  and  produces 
as  its  output  binarized  words.  This  is  to  say  that  it 
extracts  one  word  from  a  page  of  text  and  outputs  it  as  a 
black  character  on  a  white  background.  (Note  that  a  word  may 
be  one  or  more  concatenated  letters.) 

This  level  consists  of  two  algorithms,  each  of  which 
performs  a  separate  operation  on  its  input.  The  first 
algorithm  performs  a  noise  reduction  ,  ration.  It  takes  the 
input  greyscale  text  and  threshold-averages  the  levels  to 
produce  binary  text. 
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Figure  *.l  PACKED  VIDEO  FORMAT 

As  a  final  note  on  the  digitization  stage,  it  is 
important  to  remember  that  each  greylevel  in  a  video  file  is 
stored  as  an  integer  value  of  four  bits  in  a  computer  word. 
(This  becomes  an  important  fact  in  later  processing 


procedures. ) 


digitized  greyscale  pictures.  The  next  stage  would  have  to 
work  with  the  given  output. 

4.2  SYSTEM  ASPECTS 

The  digitization  equipment  consisted  of  a  light  table 
with  a  video  camera  mounted  onto  it.  Input  text  placed  onto 
the  device  is  digitized  into  a  256x256  video  file  with  16 
greylevels  (0  corresponding  to  black  and  15  to  white). 

In  order  to  understand  the  source  code  of  certain 

programs  which  follow,  it  is  important  to  mention  that  the 
digitizer  "packs"  four  pixels  into  one  computer  word  (see 

Figure  4.1).  One  computer  word  is  16  bits  long.  Thus, 

processing  requires  the  " UNPACK ”-ing  and  the  " REPACK "-ing  of 
pixel  points. 

Even  though  a  complete  digitized  picture  is  256x256 
pixels  wide,  the  actual  useable  video  size  is  only  2*0x256 
pixels  wide.  The  last  sixteen  rows  of  each  video  are  zero 
filled  to  produce  a  black  strip  at  the  bottom  of  each 
picture.  (The  reason  for  zero  filling  the  last  16  lines  is 
is  that  the  output  display  monitor  is  only  2*0x256  pixels.) 
All  65,536  pixels  of  a  video  file,  however,  are  accessible 
for  internal  processing. 

For  most  of  the  programs,  the  actual  manipulation  of 
pixel  points  was  handled  via  PICBUF.  This  is  a  library  of 
Fortran-callable  subroutines  which  hides  the  system- 
dependent  details  of  unpacking  and  repacking.  (For  more 

details  on  PICBIJF  see  Appendix  A). 


VI.  THE  RECOGNITION  STAGE 


The  following  chapter  discusses  the  software 

implementation  of  the  final  level  of  the  proposed 
recognition  algorithm.  Since  the  Recognition  stage  is  the 
actual  level  where  input  characters  are  identified,  it 
becomes  the  most  important  level  of  the  algorithm.  Enough 
code  was  generated  to  allow  for  testing  of  the  scheme.  Thus, 
in  most  cases,  the  interface  between  programs  were  not 
written.  They  would  only  allow  for  total  automatic 
processing  and  could  be  written  at  a  later  date. 

Chapter  3  described  the  overall  processing  procedure 
for  recognition.  This  chapter  will  describe  the  actual 
software  implementation.  Recall  that  identification  comes 
from  r e - c o r r e 1  a t i o n  of  energy  equalized  characters.  The 
first  section  will  discuss  the  software  for  correlation  and 
the  second  section  will  discuss  the  soft',  .'are  for  energy 
modification. 

6.1  software  developed 

The  Recognition  stage  accepts  as  input  words  produced 
by  the  Pre-processing  stage.  For  the  reasons  discussed  in 
Chapter  2,  it  performs  a  correlation  operation  on  the  word 
using  Fourier  transforms.  Chapter  3  briefly  described  the 
details  of  the  correlation  process.  Figure  6.1  shows  a 
flowchart  of  the  programs  which  implement  the  procedure. 
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Figure  6.1  FLOWCHART  OF  RECOGNITION  PROGRAMS 
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Sections  6.2  and  6.3  will  describe  in  detail  the  operation 
of  each  of  those  programs. 

As  envisioned,  an  input  word  is  first  correlated 
with  each  of  the  templates  in  the  template  set.  A  record  is 
then  kept  of  the  positions  where  each  template  peaked  with 
the  input.  At  each  peak  location,  a  box  is  centered  on  the 
peak  and  a  second  series  of  correlations  is  performed. 
During  this  second  pass,  the  energy  of  the  box  and  the 
respective  template  are  set  equal  to  one  another.  This  time, 
correlation  will  produce  peaks  of  varying  values.  The 
highest  peak  value  determines  which  template  correlated  the 
best. 

Mote  that  at  this  time  of  research,  some  of  the 
functions  just  described  are  not  fully  implemented  in 
software.  Instead  of  being  totally  automatic,  the  presented 
functions  require  that  the  user  pass  the  values  obtained 
from  one  program  to  that  of  the  next. 

6.2  ?j*22?SSS 

The  first  step  in  the  correlation  process  is  to  compute 
the  Fourier  transform  of  the  input  word.  This  is 
accomplished  by  using  the  programs  VTC  and  DIRECT.  The 
program  VTC.FR  (Video  To  Complex)  converts  a  video  file  of 
integer  numbers  into  a  complex  file  with  the  imaginary  part 
set  equal  to  zero. 

The  program  DIRECT. SV  performs  the  actual  computation 
of  the  Fourier  Transform.  It  was  a  supplied  program  which 


could  only  be  executed.  The  original  source  code  was  not 
available. 

Since  the  characters  of  the  template  set  are  known 
beforehand,  the  Fourier  transform  of  each  template  is 
computed  before  the  beginning  of  the  process.  This  computa¬ 
tion  is  performed  only  once  and  the  results  are  stored  in 
memory.  The  computation  of  the  transform  for  the  templates 
requires  an  additional  program,  QSHIFT.FR  ( Q u a r t e r S H I  FT )  .  In 
order  for  the  correlation  peaks  to  be  centered  properly  on  a 
letter,  the  template  must  be  quartered  and  the  two  diagonal 

corners  must  be  interchanged. 

The  program  which  performs  the  actual  correlation 
process  is  CMULT1.FR  (Complex  Multiplication).  From  Digital 

Signal  Processing,  recall  that  correlation  in  the  spatial 

domain  corresponds  to  multiplication  in  the  frequency 
domain.  The  program  multiplies  the  complex  value  of  a 
pixel  from  one  file  with  the  complex  conjugate  value  of 

another  file.  It  therefore,  performs  the  correlation  process 

by  complex  conjugate  multiplication  on  a  pixel  by  pixel 

level  . 

The  process  is  completed  by  inverse  Fourier 
transforming  the  correlated  file  and  converting  the  complex 

file  into  a  video  file.  The  program  INVERSE. SV  calculates 

the  inverse  transform.  The  program  CTV.fr  (Complex  To  Video) 
performs  the  conversion. 

The  program  CTV.FR  linearly  scales  the  values  of  the 


complex  file  so  that  the  highest  value  in  the  file  has  a 
value  of  15.  Therefore,  this  program  will  always  produce  a 
peak  correlation  value  of  15.  (After  the  inverse  transform 
operation,  the  video  file  contains  complex  numbers  with  the 
imaginary  part  equal  to  0.) 

After  correlation,  the  program  FPEAK.FR  (Find  Peak) 
will  determine  the  highest  correlation  peak  in  the 
correlated  scene  and  determine  the  location,  or  loactions, 
of  the  peak.  (These  values  can  be  passed  on  to  the  energy 
modification  programs.  At  this  point  in  time,  these  values 
must  be  manually  entered.) 

6 . 3  energy  modification  PROCESS 

As  discussed  in  Chapter  2,  a  correlation  peak  is 
maximum  if  the  two  shapes  have  the  same  energy  and  the  same 
shape.  After  the  first  correlation  process  determines  the 
locations  of  peaks,  the  area  centered  on  the  peak  must  be 
energy  modified  to  equal  that  of  a  template.  For  the  present 
implementation,  if  the  energy  in  the  area  is  greater  than 
that  of  a  template,  then  the  energy  of  the  area  is  set  equal 
to  that  of  the  template.  If  the  template  has  higher  energy, 
then  the  template  energy  is  set  to  that  of  the  area.  The 
program  which  performs  this  operation  is  MEND.FR  (Modify 
ENergy  Distribution). 

For  greyscale  characters,  the  energy  in  a  word  is 
computed  by  taking  the  square  root  of  the  sums  of  the 
squares  of  the  individual  pixel  values  within  a  given  area: 

49 


Since  all  of  the  characters  in  this  level  have  been 
Pre-Processed,  the  energy  becomes  the  square  root  of  the 
number  of  black  pixels  in  a  given  area. 

After  the  energies  are  established,  both  template  and 
input  are  r e - c o r r e 1  a t e d .  The  same  procedure  follows  as  in 
the  first  correlation  pass,  except  that  at  the  end,  program 
CTV5.FR  is  used.  Like  CTV.FR,  this  program  performs  a  linear 
scaling  of  the  complex  pixel  values,  but  at  a  fixed  range. 
Therefore,  each  correlation  will  produce  a  different 


correlation  peak. 


V  I  I 


RESULTS 


The  following  chapter  presents  the  results  of  some 
sample  tests  conducted  to  evaluate  the  performance  and 
limitations  of  the  proposed  recognition  scheme.  It  is 
expected  that  research  will  continue  until  a  fully 
operational  device  is  developed. 

The  first  series  of  tests  were  performed  with  letters 
from  the  same  font  style.  The  second  series  of  tests 
were  performed  with  two  differing  font  styles. 

Throughout  testing,  it  was  assumed  that  the  maximum 
width  and  height  of  a  given  character  was  known.  This 
provided  the  dimensions  for  the  window  in  the  energy 
modification  process. 

7 • 1  §ame  FONT  STYLE 

Figure  7.1  shows  the  characters  of  the  template  set 
used  for  testing.  Single  letters  and  concatenated  letters 
were  generated  from  this  template  set  and  used  for  the  first 
series  of  tests.  (Appendix  A  contains  the  source  listing  of 
the  programs  used  to  generate  the  concatenated  letters.) 

The  first  test  examined  the  recognition  scheme's 
ability  to  recognize  single  letters.  The  second  series  of 
tests  examined  its  ability  to  recognize  letters  within  a 
string  of  letters. 
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7  .  1  A  SINGLE  LETTER  TEST 

Figure  7.2  shows  a  capital  letter  F.  This  character  was 
input  into  the  Recognition  stage  for  identification.  This 
produced  a  series  of  correlation  peaks.  Figures  7.3  through 
7.7  show  the  correlation  peaks  for  the  second  pass  with  the 
letters  B,  C,  E,  F,  and  H.  The  peak  value  for  the  F 
template  was  the  highest. 

7. IB  CONCATENATED  LETTERS 

Figure  7.8  shows  five  concatenated  letters  as  the 
output  of  -he  second  stage.  Figures  7.9  through  7.13 
shows  the  results  after  the  first  correlation  pass  with  the 
"A",  "B",  "C",  "D",  and  "E"  templates,  respectively.  The 
highest  correlation  peak  at  the  first  position  in  the  word 
was  for  the  "A"  template.  The  highest  peak  at  the  second 
position  was  for  the  "B"  template.  The  highest  correlation 
peaks  for  the  third,  fourth,  and  fifth  positions  were  for 
the  "C",  "D",  and  "E"  templates,  respectively. 

Correlating  the  word  with  the  "F"  template  also 
produced  a  peak  at  the  fifth  positon  (see  Figure  7.14).  A 
window  was  therefore  centered  at  the  location.  Figure  7.15 
shows  the  peaks  for  re-correlating  the  "E"  template  with  the 
window  energy  set  equal  to  that  of  the  "E"  template.  Figure 
7.16  shows  the  same  peaks  for  the  "F"  template.  The  value  of 
the  second  correlation  peak  was  highest  for  the  "E" 
template. 
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FIGURE  7.3  INPUT  CORRELATED  WITH  TEMPLATE  B 


FIGURE  7.5 


INPUT  CORRELATED  WITH  TEMPLATE  E 


FIGURE  7.6  INPUT  CORRELATED  WITH  TEMPLATE  F 
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FIGURE  7.18  TEST  SET  CORRELATED  WITH  TEMPLATE  B 
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GHIJK 

MNOPQR 

STUVWX 


FIGURE  7.17 


TEST  SET  NUMBER  2 


7  .  2 


LARGER  FONT  STYLE 


Figure  7.17  shows  the  characters  of  another  font  style. 
These  letters  are  approximately  3  0'/.  larger  (total  amount  of 
pixels  contained  in  a  character)  than  those  of  the  first 
template  set.  This  second  template  set  was  used  as  input  to 
the  Recognition  stage. 

The  first  pass  correlation  of  the  template  "B"  (from 
the  first  template  set)  with  the  second  template  set 
produced  the  peaks  shown  in  Figure  7.18.  The  template  "B" 
correlated  well  with  every  input  character. 

From  the  discussion  in  Chapter  2,  it  appeared  as  if  the 
recognition  algorithm  might  be  helped  if  the  target  letters 
were  not  so  thick.  Figure  7.19  shows  a  thinned  version  of 
the  second  template  set.  (The  program  VTHIN.FR  was  used  to 
thin  the  letters.  The  source  code  listing  for  this  program 
is  in  Appendix  A.)  The  template  " B "  of  the  first  set  was 
then  correlated  with  the  thinned  version.  Figure  7.20  shows 
the  resulting  peaks.  The  template  "B"  produced  a  peak  with 
the  letter  G. 

Since  the  first  correlation  pass  did  not  produce  a 
peak  at  the  proper  loacation  testing  was  not  continued.  It 
appears  as  if  the  "template  height  matching"  function  of  the 
algorithm  must  be  implemented.  However,  time  limitations  did 
not  allow  for  this. 
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FIGURE  7.15  CORRELATION  WITH  ENERGY  SET  TO  TEMPLATE  E 
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CORRELATION  WITH  TEMPLATE  C 
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VIII.  DISCUSSION 


Those  Casks  which  seem  trivial  to  a  human  seem  to  be 
impossible  for  a  computer.  With  all  their  computing  power 
and  megabytes  of  memory,  the  digital  computer  is  in  reality 
a  very  stupid  machine.  It  can  only  perform  those  tasks  which 
a  person  allows  it  to  perform  (by  writing  the  proper 
program) . 

8 • 1  CONCLUSION 

This  thesis  examined  the  use  of  a  two  pass  correlation 
technique  for  character  recognition.  The  first  pass  would 
produce  a  series  of  correlations  which  showed  which 
templates  matched  best  with  the  input.  During  the  second 
pass,  the  energy  of  both  input  and  template  would  be  set 
equal  to  one  another.  Correlation  would  then  produce  a 
maximum  peak  for  inputs  with  the  same  shape. 

Testing  with  a  highly  constrained  font  set  showed  that 
the  algorithm  would  work  for  isolated  letters  as  well  as  for 
concatenated  letters  if  the  template  set  was  derived  from 
the  same  font  style  as  that  tested. 

Testing  with  a  template  set  different  from  that  of  the 
input  test  set  showed  that  the  algorithm  was  very  sensitive 
to  size  variation.  Thinning  letters  to  match  the  template 
pixel  thickness  seemed  to  aid  the  process.  However, 
recognition  was  still  not  possible.  It  appears  as  if  the 


height  of  the  template  must  also  be  adjusted.  Thus,  energy, 
thickness,  and  height  must  all  be  considered. 

8.2  RECOMMENDATIONS 

It  is  apparent  that  much  more  research  must  be 
performed  before  a  recognition  is  developed  which  is  truely 
font  invariant.  For  the  algorithm  presented,  three 
recommendations  are  made  for  future  research: 


1.  Implement  in  software  the  function  which  adjusts  the 


height 

of  a 

template 

t  o 

match  that  of  the  input. 

The 

program 

LINES 

could  be 

e  a  s 

ily  modified  to  calculate 

the 

average 

he i g h  t 

of  a  line 

o  f 

print. 

2.  Testing  should  continue,  but  for  a  greater  number  of 
characters.  These  tests  should  contain  both  the  upper  and 
lower  case  characters.  (Since  the  height  was  not  adjusted  in 
this  thesis,  capital  letters  did  not  correlate  well  with 
lower  case  letters.) 

3.  Try  correlation  with  the  Fourier  transform  clipped  after 
the  third  harmonic.  Most  of  the  information  content  in  a 
scene  is  contained  in  the  low  order  harmonics.  Serifs  and 
jagged  edges  contribute  principlely  to  the  high  order 


ha  r  mon i c  s  . 
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This  part  of  the  thesis  contains  the  software  developed 


for  the  pre-processing  and  recognition  stages.  Each  program 
is  listed  in  alphabetical  order  with  the  exception  of  DIRECT 
and  INVERSE.  They  were  supplied  as  executable  programs 
without  source  code. 

From  documentation  in  the  Signal  Processing  Lab,  DIRECT 
is  a  Fortran  program  written  by  R.  W.  Schafer  on  30  June  78. 
The  purpose  of  the  program  is  to  compute  a  2-D  forward  FFT 
on  disk  using  the  Eklundh  algorithm.  The  maximum  input  array 
size  is  312x512.  Input  must  be  in  binary  format  and  contain 
complex  data  with  the  imaginary  part  equal  to  zero. 

INVERSE  is  another  Fortran  program  written  by  Schafer 
on  15  Feb  77.  This  program  computed  the  2-D  inverse  FFT  on 
disk  using  the  Eklundh  algorithm.  The  inverse  transformed 
array  is  normalized  by  division  with  N**2. 

Two  notes  are  required  before  this  thesis  concludes. 
The  first  is  that  the  program  MEND.FR  was  written  using 
integer  values  for  output  because  the  OCTEK  only  accepted 
integer  values. 

The  final  note  is  that  line  thinning  was  examined  in 
the  later  part  of  this  thesis.  The  program  VTHIN.FR  was 
written  to  perform  this  task.  The  program  RLNX.FR  must  be 
used  first  to  determine  the  maximum  number  of  pixels  which 
can  be  deleted  from  the  character. 
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PROGRAM:  CLEAN. FR 


LANGUAGE:  FORTRAN  V 


DO  4  1=1,240 

CALL  GROW( I BUF , I , 1 , 2 5 6 , I ARRA Y ) 
DO  6  J=  1  ,  256 

K= I  ARRAY ( J )  +  2 
AHIST(K)=AHIST(K)+1 . 0D0 
CONTINUE 
CONTINUE 


GET  EACH  PIXEL 
VALUE 
(240  ROWS, 

256  COLUMNS) 
CALCULATE  THE 
HI STOGRAM 


DO  8  1=1,18 

I F ( AH  I ST(  I  )  . EQ . 0 ) K=K+ 1 
CONTINUE 

I  F ( K . NE .  17  )  GOTO  8  9 

TYPE "THERE  IS  ONLY  ONE  LEVEL  OF  BRIGHTNESS" 
GOTO  99 

I  F ( K  .  NE .  16  )  GOTO  90 

TYPE "THERE  ARE  ONLY  TWO  LEVELS  OF  BRIGHTNESS" 
GOTO  99 


DETERMINE 
IF  VIDEO 


BLANK 


B INARY 


C  =  1 

DO  10  1=2,1 


IF(  ( AH  I ST(  I  )  . GT . AH  I ST(  1-1  )  ) 

.  AND  . 

(AHIST(I) .GE.AHIST(I+1)))G0T0  91 


GOTO  10 

AMAX(C)=AHIST( I ) 
C  =  C+  1 


DETERMINE 

THE 

PEAKS 


CONTINUE 


DO  451  1=1,18 

WRITE  FREE ( 10)"AHIST",I ,AHIST( I) 
1  CONTINUE 


80 


c 

c 


D  =  C  -  2 


DO  12  I  -  1  ,  D 
J  =  C  -  I 

DO  14  K  = 1  , J 

I F ( AMAX ( K ) . GE . AMAX ( K+l ) )GOTO 
T  =  AM A  X ( K ) 

AMAX ( K ) =  A  M A  X ( K  + 1  ) 

AMAX ( K+l  ) =  T 
CONTINUE 
CONTINUE 


BUBBLE 


PEAK 

VALUES 


DO  16  1=2,17 

I F ( AH  I ST(  I  )  . NE . AMAX (  1  )  )GOTO  16 
AN  DX  =  I 
GOTO  991 
CONTINUE 


DETERMINE 
GREYLEVEL 
CORRESPONDING 
TO  THE  FIRST 
PEAK  VALUE 


N  =  D+  1 
J  =  2 

DO  18  1=2,17 

I  F ( AH  I ST(  I)  .  NE . AMAX( J )  ) GOTO  18 
N  =  N  -  1 
B  N  D  X  =  I 

IF(  (  I  .  GT .  7 )  -  AND-  (N.GT.O)  )GOTO  92 
GOTO  992 
J  =  J+l 
CONTINUE 


DETERMINE 
GREYLEVEL 
CORRESPONDING 
TO  THE  SECOND 
PEAK  VALUE 


I F ( BNDX . LT . ANDX ) GOTO  993 
E  =  A  N  D  X 
ANDX=BNDX 
B  N  D  X  =  E 


ENSURE  THAT 
ANDX  >  BNDX 


C  =  1 

DO  20  I=BNDX,ANDX 


IF((AHIST(I).LE.AHIST(I-1)) 

.  AND  . 

( AH  I ST(  I  )  .  LT . AH  I  ST<  1  +  1  )  )  )GOTO  93 


GOTO  20 

AM  I N ( C ) = AH  I  ST (  I  ) 
C  =  C+  1 


CONTINUE 


DETERMINE 


VALLEYS 


C 


ro  ro 


c 


c 


I F ( C . EQ . 2 ) GOTO  994 


DETERMINE  IF  C 

THERE  ARE  ONLY  C 

ONLY  TWO  C 

VALLEYS  C 


D  =  C  -  2 

DO  2  2  I  =  L  ,  D 
J  =  C  -  I 


DO  2  4  K  =  t  ,  J 

I F ( AMIN ( K ) . LE . AMI N ( K+ 1 ) ) GOTO  24 
T  =  AM I N ( K ) 

AMIN ( K ) =AMI N ( K  + 1 ) 

AM  I N ( K+  1  ) =T 
CONTINUE 


CONT  I  NL'E 


BUBBLE 


VALLEYS 


DO  26  1=2,17 

IF(AHIST(I) .NE.AMIN( I ))GOTO 
MN  DX  =  I 
GOTO  995 
CONTINUE 


DETERMINE 
GREYLEVEL 
CORRESPONDING 
TO  THE  FIRST 
VALLEY 


MNDX=MNDX-2 

IF(MNDX.GT.  i  U  )  HffnUUKt  CANNOT 
I F ( MNDX . GT . 1 0 ) GOTO  9999 


BE  CLEANED' 


DETERM  I  NEC 
IF  C 

VIDEO  C 

CAN  BE  C 

PROCESSEDC 


M  I  D  =  I  NT (  (ANDX  +  BNDX)/2) 


CALCULATE 
MIDPOINT 
BETWEEN  THE 
TWO  PEAKS 


WRITE  FREE(  l  0  )  "ANDX-"  , ANDX-2 
WRITE  FREE (  I  0  )  " BNDX  =  ”  ,  BNDX- 2 
WRITE  FREE (  1  0  )  "MNDX  =  ", MNDX 
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DO  55  I RO W= 
CALL 


,  240 

GROW(  IBUF,  I  ROW 


256 , IARRAY) 


GET  A  ROW 
OF  PIXELS 


NT  =  0 

DO  60  I =  MN  D  X  ,  10 
NT  =  NT+ 1 

DO  65  J  =  l  ,256 

I F ( I  ARRAY ( J  )  . LE .  I  )GOTO  95 
IARRAY(J)=IARRAY(J)+IWHT 
GOTO  65 

IARRAY(J)=IARRAY(J)+IBLK 

CONTINUE 

CONTINUE 


THRESHOLD 
POINT  VARIES 
BETWEEN  THE 
FIRST  VALLEY 
AND  10 

CLIP  PIXELS 
BELOW 
THRESHOLD 
TO  0  AND 
SUM  PIXEL 
VALUES 


DO  70  J  =  1  ,256 

IARRAY ( J ) = I NT( IARRAY ( J) /NT) 
I F (  I  ARRAY ( J )  . LE . MI D )GOTO  96 
I  ARRAY ( J )  =  IWHT 
GOTO  70 

IARRAY(J)=IBLK 

CONTINUE 


AVERAGE 

THRESHOLDED 

VALUE 

CLIP  VALUES 
BELOW 
MIDPOINT 
TO  0 


CALL  PROW (IBUF,  IROW,  1  ,256,  IARRAY) 


CONTINUE 


DO  75  1=1,256 

IARRAY ( I ) = 1 5 
CONTINUE 

DO  85  1=241,256 

CALL  PROW(  I BUF ,  I  ,  1  ,  2  5  6  ,  I  ARRAY  ) 
CONTINUE 


BLANK  OUT 
LAST  16  ROWS 
OF  THE  VIDEO 
FILE 


CALL  PICOUTC IBUF, INAME, 1 ) 


TYPE" < 7 >F I N I  SHED  :  PLEASE  CALL  SPOTS" 


PICBUF  OUTPUT 
BLOCK 


CALL  RELB(IBUF) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


STOP 


Q***************************************************************£ 

C  C 

c  c 

C  ORIGINAL  PROGRAM:  CPTOVD.FR  BY  R.  WELLS  C 

C  C 

C  PROGRAM:  CTV.FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  PROGRAM  FINDS  THE  MAXIMUM  OF  A  256X256  C 

C  COMPLEX  FILE  AND  GIVES  THE  LOCATION  OF  THE  C 

C  MAXIMUM.  ALSO,  THE  PROGRAM  CONVERTS  THE  FILE  C 

C  TO  PACKED  VIDEO  FORMAT  BY  PERFORMING  A  LINEAR  C 

C  SCALING  TO  THE  RANGE  OF  0  -  15  AND  TRUNCATING.  C 

C  C 

C  C 

C  C 

£***********************************************:****************q 

c  c 

c  c 


INTEGER  I FLNM( 7 ) , OFLNM( 7 ) 

INTEGER  TEMP (  5  12)  ,0UT(2048)  ,  IMAX , JMAX , MAXCOL , MAXROW 
COMPLEX  IN( 2048) 

REAL  RMAG . MAXVAL , MI NVAL , RANGE 


C  ;  C 

ACCEPT  "ENTER  INPUT  FILENAME:  "  ;  C 

READ  (11,1995)  IFLNM(l)  ;  ENTER  C 

1995  FORMAT  (S40)  ;  INPUT  AND  OUTPUT  C 

;  FILENAMES  C 

ACCEPT  "ENTER  OUTPUT  FILENAME:  "  i  C 

READ  (11,1995)  OFLNM(l)  ;  C 

C . ; . C 

C  ;  C 

OPEN  1 , I FLNM  ;  C 

;  C 

CALL  DFILW  (OFLNM.IER)  ;  OPEN  C 

IF( lER.Nfc. 1 )GOTO  90  ;  INPUT  AND  OUTPUT  C 

90  CALL  CF I LW( OFLNM , 2 , KER )  ;  FILES  C 

CALL  CHECK(KER)  ;  C 

OPEN  2 , 0 F L NM , ATT - " 0 R " , E R R = 1 0 0  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

100  MAXVAL- 0.0  ;  C. 

MINVAL-0.0  ;  C 

;  C 

0031-0,31  ;  C 

;  C 

CALL  RDBLK(  1  ,  3 2  *  I  ,  IN ,  3 2  ,  I  ER  )  ;  C 

CALL  CH  E  CK (  IER)  ;  C 

;  C 

DO  3  J- 1 , 2048  ;  C 

;  DETERMINE  C 


84 


RMAG-CABS ( IN ( J  )  ) 

I F ( RMAG . LE . MAXVAL)GO  TO  3 
MAXVAL-RMAG 
I  M  A  X  =  I 
JMAX»J 


3  CONTINUE 


MAXVAL 

AND 

MINVAL 


RANGE* MAXVAL-MINVAL 

MAXCOL- MO D(J MAX-1 ,256)+! 
MAXROW»8*IMAX+INT(  (JMAX-1  )/256)+l 


CALCULATE 

RANGE 

MAXCOL 

MAXROW 


TYPE" MAX 

COL 

XX 

M 

,  MAXCOL 

TYPE "MAX 

ROW 

3 

II 

,  MAXROW 

TYPE "MAX 

VAL 

a 

ft 

,  MAXVAL 

TYPE "MIN 

VAL 

a 

II 

.MINVAL 

DO  5  1*0,31 

CALL  RD8LK(  1  ,  32* I  ,  IN  ,  32  ,  IER) 

CALL  CHECK ( I ER ) 

DO  4  J  *  1  ,2048 

OUT(J)-ANINT(15.0*(CABS(IN(J)) 
/  -MINVAL) /RANGE) 

4  CONTINUE 

CALL  REPACK( 5 1 2 , OUT , TEMP ) 

CALL  WRBLK ( 2  ,  2  *  I  , TEMP , 2  ,  I ER ) 

CALL  CHECK(IER) 

5  CONTINUE 


SCALE  VALUES 
TO 

0-15 


CALL  RESET 
STOP"  <7X7: 


E 


c  c 
c  c 
c  c 

C  ORIGINAL  PROGRAM:  CPTOVD.FR  BY  R.  WELLS  C 
C  C 
C  PROGRAM:  CTV5.FR  LANGUAGE:  FORTRAN  V  C 
C  C 
C  WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE  C 
C  C 
C  PURPOSE : TH I S  PROGRAM  FINDS  THE  MAXIMUM  AND  MINIMUM  OF  A  C 
C  256X256  COMPLEX  FILE  AND  GIVES  THE  LOCATION  OF  THE  C 
C  MAXIMUM.  THE  PROGRAM  CONVERTS  THE  FILE  TO  PACKED  C 
C  VIDEO  FORMAT  BY  PERFORMING  A  LINEAR  SCALING  OF  C 
C  0.0  -  50000.0  TO  RANGE  OF  0  -  15.  C 
C  C 
C  C 
C  C 

Q***************************************************************^ 

c  c 
c  c 


INTEGER  I FLNM( 7  )  , OFLNM( 7  ) 

INTEGER  TEMP(5I2) ,0UT(2048) , IMAX , JMAX , MAXCOL , MAXROW 
COMPLEX  IN(2048) 

REAL  RMAG , MAXVAL , MINVAL , RANGE 


C  ;  C 

ACCEPT  "ENTER  INPUT  FILENAME:  "  ;  C 

READ  (11,1995)  IFLNM(l)  ;  ENTER  C 

1995  FORMAT  (S40)  ;  INPUT  AND  OUTPUT  C 

;  FILENAMES  C 

ACCEPT  "ENTER  OUTPUT  FILENAME:  "  ;  C 

READ  (11,1995)  OFLNM( 1 )  ;  C 

C 
C 

C  ;  C 

MAX VAL« 5 0000 . 0  ;  SET  MAXVAL  C 

C  ;  C 

C . ; . C 

C  ;  C 

OPEN  1,1 FLNM  ;  C 

;  C 

CALL  DFILW  (OFLNM.IER)  ;  OPEN  C 

I F ( I ER . NE . 1 ) GOTO  90  ;  INPUT  AND  OUTPUT  C 

90  CALL  CF I LW( OFLNM ,  2  ,  KER  )  ;  FILES  C 

CALL  CHECK(KER)  ;  C 

OPEN  2  , OFLNM , ATT*"OR"  ,  ERR- 1 00  ;  C 

C  ;  C 

C . ; . C 
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n  n 
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F/G  5/7 


UNCLASSIFIED 


NL 


C . ; . C 

C  ;  C 

C  ;  C 

100  MINVAL-0.0  ;  INITIALIZE  C 

C  ;  MINVAL  C 

C  ;  C 

C . ; . C 

C  ;  C 

RANGE* 50000 . 0  ;  INITIALIZE  C 

C  ;  RANGE  C 

C . ; . C 

C  ;  C 

TYPE "MAX  VAL  *  "  , M A  X  V  A  L  ;  C 

TYPE "MIN  VAL  -  ".MINVAL  ;  I/O  C 

C  ;  C 

C . ; . C 

C  ;  C 

DO  5  I =0 , 3 1  ;  C 

;  C 

CALL  RDBLK(  1  ,  32*1  ,  IN  ,  32  ,  IER)  ;  C 

CALL  C  H  E  CK (  IER)  ;  C 

;  C 

DO  4  J " 1 ,2048  ;  SCALE  VALUES  C 

OUT ( J ) » AN INT(  1 5 . 0* ( CABS ( I N ( J  )  )  ;  C 

/  -MINVAL) /RANGE)  ;  TO  C 

4  CONTINUE  ;  0-15  C 

;  C 

CALL  REPACK( 5 12 .OUT, TEMP)  ;  C 

CALL  WRBLK( 2  , 2* I  , TEMP  ,  2 ,  IER ) 

CALL  CH  E  C  K (  IER) 

;  C 

5  CONTINUE  ;  C 

C  ;  C 

C . ; . C 

c  c 

CALL  RESET 
END 
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O  O 


*********************************************£ 

C 

c 

PROGRAM:  FCHAR.FR  LANGUAGE:  FORTRAN  V  C 

C 

WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE  C 

C 

PURPOSE:  THIS  PROGRAM  WILL  SEPARATE  EACH  CHARACTER  FROM  C 

THE  NEXT  IF  THERE  EXISTS  A  SPACE  BETWEEN  THEM.  C 

C 

************************★***********************************£ 

c 


IMPLICIT  INTEGER(C , D , E) 

INTEGER  IARRAY(512) ,  I BUF ( 3  00  ) ,INAME(20) 
INTEGER  DCNTR(20),DNTR(20),ECNTR(20),ENTR(20) 
INTEGER  I  C S TK (  100) 


ACCEPT" ENTER  FILENAME:  " 
READ (  1  1  ,  199)1  NAME (  1  ) 
FORMAT ( S  4  0  ) 


INPUT 
F I LENAME 


CALL  P ICFMT( I BUF ,  INAME ,  I  FLG )  ; 

CALL  GFMT( I BUF , NROWS , NCOLS , NBITS , IMODE )  ;  PICBUF  INPUT 

;  BLOCK 


CALL  MAKB ( I BUF ) 

CALL  P ICIN( I BUF , INAME , IFLG ) 


OPEN  1 , " LD AT . TX" 

READ  B I  NARY (  l  ) DCNTR , ECNTR , C , LIN 
OPEN  2 , "WDAT.TX" 

READ  B INARY ( 2 ) DNTR , ENTR , D , I WORD 


READ  LINE 
AND  WORD 
POSITIONS 


TYPE "LINE:", LIN,"  WORD:", I  WORD 


I/O 


TYPE"FCHAR  WORKING" 


I S  C  0  L  *  D  NTR ( IWORD) 
ISCOL»ISCOL+l 


CALCULATE 
STARTING  COLUMN 


CALL 


P I C I N (  IBUF,  INAME,  IFLG) 


GET 


VIDEO 


FILE 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


o  o 


c . 

COL* I SCOL 

E-l  ; INITIALIZE  COUNTERS  C 

ICSTK(E)=COL  ;  TO  STARTING  C 

;  POSTIONS  C 

DC-DCNTR( LIN )  ;  C 

EC-ECNTR(LIN)  ;  C 

C . ; . C 

C  ;  C 

TYPE" START  COL-  ",ICSTK(E)  ;  C 

;  I/O  C 

TYPE" START  ROW-  ",DC-1  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

DC-DC-l  ;  C 

IF ( DC . LT .  1  )  DC- 1  ;  CALCULATE  START  C 

EC-EC+L  ;  AND  END  OF  ROW  C 

IF( EC . GT . 2 5 6 ) EC- 2 5 6  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

IST-DC  ;  C 

IEN-EC  ; DETERM I N  E  IF  PIXEL  C 

1  DO  10  I  ROW- DC , EC  ;  IS  BLACK  C 

CALL  GPNT(  I  BUF  ,  IROW , COL ,  I  VAL  )  ;  C 

IF( I VAL . EO . 0 )GOTO  100  ; IF  WHITE-->  C 

ICSTK(E)-COL  -.CONTINUE  DOWN  COLUMN  C 

E-E+l  ;  C 

10  CONTINUE  ; I F  BLACK-->  C 

; FIND  ANOTHER  PATH  C 

GOTO  500  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

100  DC  - IROW  ;  C 

ITEMP-COL  ;  C 

I R - I  ROW  ;  FIND  CLOSEST  C 

IFR-ITEMP  ;  BLACK  PIXEL  TO  C 

200  IFR-IFR+1  ;  TO  THE  RIGHT  C 

CALL  GPNT (  I BUF  ,  I  R  ,  IFR , I VAL )  ;  C 

I F ( I VAL . EQ . 0 ) GOTO  200  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

C  ;  C 

IB-ITEMP  ;  FIND  CLOSEST  C 

300  IB-IB-l  ;  BLACK  PIXEL  TO  C 

CALL  GPNT( I BUF , I R , I B , I VAL )  ;  TO  THE  LEFT  C 

I F ( I VAL . EQ . 0 )GOTO  300  ;  C 

C  ;  C 

C . . . C 


o  n 


C . ; . C 

C  ;  C 

MR- I R-  l  ;  C 

DO  964  I TT - 1 B ,  I  TEMP  ;  DETERMINE  IF  LEFT  C 

CALL  GPNT(  I BUF , MR ,  ITT  ,  I VAL  )  ;  BLACK  PIXEL  IS  C 

I F ( I VAL . EQ . 0 )GOTO  186  ;  CONNECTED  TO  C 

964  CONTINUE  ;  PREVIOUS  PIXEL  C 

C  ;  C 

C . ; . C 

C  ;  C 

IFD-IABS( IFR-ITEMP)  ;  C 

I BD- I ABS ( ITEMP- I B )  ;  DETERMINE  C 

COL- I F  R  ;  FURTHEST  BLACK  C 

I F( IFD . GT .  I BD )COL=  I  B  ;  PIXEL  AND  NEW  C 

I CSTK ( E ) -CO  L  ;  STARTING  COLUMN  C 

C  ;  C 

C . ; . C 

C  ;  C 

GOTO  1  ;  REPEAT  PROCESS  C 

C  ;  WITH  NEW  COLUMN  C 

C . ; . C 

C  ;  C 

500  DO  8  1-1,256  ;  C 

I  ARRAY (  I  ) » 1  5  ;  C 

8  CONTINUE  ;  BLANK  OUT  C 

OTHER  C 

J-l  ;  CHARACTERS  C 

DO  20  I-IST.IEN  ;  TO  RIGHT  C 

ILEN«256-ICSTK( J)  +  l  ;  C 

CALL  PROW( I BUF ,  I  ,  ICSTK( J )  ,  I LEN ,  I  ARRAY )  ;  C 

J  -  J+- 1  ;  C 

20  CONTINUE  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

CTOT-O  ;  C 

I LEN- I ABS ( ENTR (  I  WORD ) -ONTR ( I  WORD )  )  ;  C 

DO  30  I»IST;IEN  ;  DETERMINE  IF  C 

CALL  GROW(  I BUF ,  I  , DNTR <  IWORD )  ,  I LEN  ,  I  ARRAY )  ;  PROCESS  HAS  C 

DO  40  J-l.ILEN-l  ; BLANKED  ENTIREC 

I F ( I ARRA Y ( J ) . EQ . 0 ) CTOT-CTOT+ l  ;  LINE  C 

4 0  CONTINUE  ;  C 

30  CONTINUE  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

C  ;  IF  LINE  IS  C 

I F ( CTOT . LT . 30 )GOTO  186  ;  BLANK  C 

C  ; REPEAT  PROCESS  C 

C . ; . C 

C  ;  C 

OPEN  3,"CDAT.TX"  ;WRITE  BLANKING  C 

WRITE  BINARY( 3) ICSTK  ;  COLUMNS  INTO  C 

C  ;  CDAT.TX  C 

C . ; . C 
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CALL  P I COUT( I BUF , INAME , I FLG ) 
TYPE"<7>FINISHED" 


PICBUF  OUTPUT 
BLOCK 
AND 
I/O 


CALL  RELB( IBUF) 


STOP 


PROGRAM:  FLINE.FR 


LANGUAGE:  FORTRAN  V 


i  v  i  •  i  ■  \i  »  5  ■  » 


**************£ 

c 

c 


c 
c 

WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE  C 

C 

PURPOSE:  THIS  PROGRAM  WILL  COPY  ONE  LINE  OF  PRINTED  TEXT  C 

(WHICH  THE  USER  HAS  CHOSEN)  AND  BLANK  OUT  THE  C 

OTHER  LINES.  C 


C 

C 

************************************************************£ 

c 

c 


IMPLICIT  I NTEGER ( C , D , E  ) 

INTEGER  IARRAY(5  12)  ,  I BUF(  300  )  ,IBUF2(300)  ,  INAME(20) 
INTEGER  DCNTR( 20) , ECNTR( 20) 


ACCEPT "ENTER  FILENAME: 
R  E  A  D (  11,199) I N  AME (  1  ) 
FORMAT ( S  4  0 ) 


OPEN  1 , "LDAT.TX" 

READ  BINARY ( 1 ) DCNTR , ECNTR , C 


ACCEPT "WHICH  LINE  TO  PROCESS:  ",LIN 
I F (  (LIN.GT.C-1  )  .OR.  (LIN  .LE.O)  )GOTO  299 


CALL 

P  I  C  F  MT (  IBUF,  INAME,  IFLG) 

CALL 

GFMT( IBUF , NROWS , NCOLS ,NBITS , 

IMODE  ) 

CALL 

MAKB ( IBUF) 

CALL 

PICIN(IBUF, INAME, IFLG) 

CALL 

PFMT(IBUF2, NROWS, NCOLS, NBITS 

, IMODE ) 

CALL 

MAKB ( I BUF2 ) 

TYPE "LINE:", LIN 
TYPE"FLINE  WORKING 


INPUT 
F I LENAME 


READ  LINE 
POSITIONS  FROM 
LD AT . TX 


INPUT 

LINE  NUMBER 


PICBUF  INPUT 
BLOCK 


I/O 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


! 


I 


I 


o  n 


. c 

;  C 

DO  51-1,  256  ;  C 

I  ARRAY (  I  )  —  1  5  ;  CREATE  A  C 

5  CONTINUE  ;  WHITE  C 

;  BCKGROUND  C 

DO  10  I- 1 , 256  ;  C 

CALL  PROW(  I  BUF2  ,  I  ,  l  ,  2  5  6  ,  I  ARRAY  )  ;  C 

10  CONTINUE  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

DO  15  I -DCNTR ( L I N  )  , ECNTR ( L IN )  ;  PUT  ORIGINAL  C 

CALL  GROW(  I  BUF  ,  I  ,  1  ,  2  5  6  ,  I  ARRAY  )  ;  LINE  ONTO  C 

CALL  PROW(  I  BUF2  ,  I  ,  l  ,  2  5  6  ,  I  ARRAY )  ;  WHITE  C 

15  CONTINUE  ;  BACKGROUND  C 

C  ;  C 

C . ; . C 

C  ;  C 

CLOSE  1  ;  C 

OPEN  1 , " LDAT . TX "  ;  WRITE  LINE  C 

WRITE  B I  NARY (  1  ) DCNTR , ECNTR , C , L I N  ;  DATA  INTO  C 

C  ;  LDAT.TX  C 

C  ;  C 

C . ; . C 

C  ;  C 

CALL  P ICOUT (  I BUF2  ,  I  NAME ,  I FLG  )  ;  PICBUF  OUTPUT  C 

CALL  RELB(IBUF)  ;  BLOCK  C 

CALL  RELB(  I BUF2  )  ;  C 

C  ;  C 

C . ; . C 

c  c 

STOP"<7>  " 

END 


£***********************************************&***************£ 

C  C 

c  c 

C  PROGRAM:  FPEAK.FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN:  15  DEC  8*  SYSTEM:  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  PROGRAM  WILL  FIND  THE  MAXIMUM  CORRELATION  C 

C  PEAK  VALUE  IN  A  VIDEO  FILE.  C 

C  C 

C  C 

£***************************************************************£ 
C  C 

c  c 

INTEGER  TEMP (  5 1 2 )  .  IN (  2048  ) 

INTEGER  I FLNM /  7  )  ,IMAX(  10)  ,JMAX(  10)  ,MAXCOL(  10)  ,MAXROW(  10) 

C  C 

c . c 

C  ;  C 

ACCEPT  "ENTER  INPUT  FILENAME:  "  ;  C 

READ  (11,1995)  IFLNM(l)  ;  ENTER  INPUT  C 

1995  FORMAT  (S40)  ;  FILENAME  C 

C . ; . C 

OPEN  l.IFLNM  ;  I/O  FILE  CONTROL  C 

C . ; . C 

C  ;  C 

MAXVAL-0  ;  C 

;  C 

DO  5  I -0 , 3  1  ;  C 

CALL  RDBLK(  1  ,  2*1  .TEMP  ,  2  ,  IER)  ;  C 

CALL  CHECK ( IER)  ;  C 

CALL  UNPACK/ 5 1 2 .TEMP , IN)  ;  DETERMINE  C 

;  MAXIMUM  C 

DO  10  J  =  1  ,2048  ;  PEAK  VALUE  C 

I F (  I  N ( J  )  .  LE . MAXVAL ) GOTO  10  ;  C 

MAXVAL-IN(J)  ;  C 

10  CONTINUE  ;  C 

;  C 

5  CONTINUE  ;  C 

C . ; . C 

C  ;  C 

K-  1  ;  C 

DO  1  5  I -0 , 3  1  ;  C 

CALL  RDBLK<  1  ,  2* I  .TEMP  ,  2  ,  I ER  )  ;  C 

CALL  CHECK/ IER)  ;  C 

CALL  UNPACK/  5 1  2  .TEMP  ,  IN)  ;  C 

;  C 

DO  20  J  »  l  ,2048  ;  DETERMINE  C 

I F ( IN ( J ) . NE . MAX VAL ) GOTO  20;  LOCATION  C 

IMAX ( K  )  -  I  ;  OF  C 

JMAX ( K  )  * J  ;  PEAK  VALUES  C 

K-K-t-1  ;  C 

20  CONTINUE  ;  C 

;  C 

15  CONTINUE  ;  C 
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» 


DO  30  I  -  1  ,  K-  1 

MAXCOL(I)*MOD(JMAX(I)-l ,256)+l 
MAXRO W (  I)«8*IMAX(  I )  +  INT( ( JMAX(  I)-l  )  /  2  5  6  )  +  1 
30  CONTINUE 
C 

C . 

C  ; 

WRITE  FREE (  1 0 ) "MAX VAL  IS:  "  , M A  X  V  A  L  ; 

T Y  P  E " AT  POSITIONS"  ; 


•♦0 
C  .  . 

c 


DO  «0  I  -  I  , K-  1 

WRITE  FREE (  1 0 )" ROW  :  ",MAXROW{  I  ) 
WRITE  FREE < 1 0 ) "COLUMN :  ”,MAXCOL(I) 


I/O 


CONTINUE 


CALL  RESET 

STOP  " <  7  >  FINISHED" 

END 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


I  c  • 

I 
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£***************************************************************£ 

c  c 

c  c 

C  PROGRAM:  FWORD.FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  PROGRAM  WILL  COPY  ONE  "WORD"  (WHICH  THE  USER  C 

C  HAS  CHOSEN)  AND  BLANK  OUT  THE  OTHER  WORDS  WITHIN  C 

C  THE  LINE.  C 


£***************************************************************£ 

C  C 

IMPLICIT  I NTEGER ( C , D , E ) 

INTEGER  IARRAY(512),IBUF(300),IBUF2(300),INAME(20) 

INTEGER  DCNTR ( 20 )  ,DNTR(20)  ,ECNTR(20)  ,ENTR(20) 


ACCEPT" ENTER  FILENAME: 
READ (  1  1  ,  199)  I  NAME (  1  ) 
FORMAT( S  4  0 ) 


INPUT 
F  I  LENAME 


OPEN  2 , " WD AT . TX" 

READ  B  I  NAR Y ( 2  ) DNTR  ,  ENTR  ,  0 


READ  WORD 
POSITIONS  FROM 
WD  A  T . TX 


ACCEPT "WHICH  WORD  TO  PROCESS:  " , IWORD 
I F (  (  IWORD . GT . D-  1  ) 

/  .  OR  . 

/  (  IWORD . LE . 0 )  )GOTO  299 


INPUT 

WORD  NUMBER 


CALL  P  I  CFMT (  I  BUF  .  I  NAME  ,  I FLG ) 

CALL  G  FMT (  I BUF , NROWS , NCOLS , NB I TS , I  MODE ) 

CALL  MAKB(IBUF) 

CALL  P  I  C I N (  I B  U  F ,  I  N  A  M  E  ,  I F  LG ) 

CALL  PFMT '  I  BUF  2  , NROWS  ,  NCO  LS  ,  N  B  I  TS  ,  I  MODE ) 

CALL  MAK  B (  I  BU  F  2  ) 


PICBUF  INPUT 
BLOCK 
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c . c 

C  ;  READ  LINE  C 

OPEN  1 , " LDAT . TX"  ;  POSITION  FROM  C 

READ  B I NAR Y (  1  ) DCNTR , ECNTR , C ,  L  I  N  ;  LDAT  .  TX  C 

C  ;  C 

C . ; . c 

C  ;  C 

C  ;  C 

TYPE"LINE:",LIN,"  WO R D  :  "  ,  I  WO R D  ;  C 

;  I/O  C 

TYPE"FWORD  WORKING"  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

DO  51=1, 256  ;  C 

I  ARRAY (  I  )  =  l  5  ;  C 

5  CONTINUE  ;  CREATE  A  C 

;  WHITE  C 

DO  10  1-1,256  ;  BACKGROUND  C 

CALL  PROW(  I BUF2  ,  I  ,  1  ,  25  6  ,  IARRAY  )  ;  C 

10  CONTINUE  ;  C 

C  ;  C 

C . ; . 

C  ;  P’ 

I LEN-ENTR (  I  WORD ) -DNTR (  I WORD )  +  1  ;ORIG  'L  C 

DO  15  I -DCNTR ( L I N ) , ECNTR ( LIN )  ;  WORD  C 

CALL  GROW(  I BUF  ,  I  ,  DNTR(  IWORD ) -  1  ,  I LEN  ,  I  ARRAY )  ;  ONTO 

CALL  PROW(  I  BUF2  ,  I  ,  DNTR (  I  WORD ) -  1  ,  I LEN ,  I  ARRAY  )  ;  WHITE  C 

15  CONTINUE  |  BACK-  C 

C  ;  GROUND  C 

C . ; . C 

0  ;  C 

CLOSE  2  ;  WRITE  WORD  C 

OPEN  2 , " WD AT . TX "  ;  POSITIONS  INTO  C 

WRITE  B  I  NARY (  2 ) DNTR , ENTR ,  D  ,  I  WORD  ;  WD AT . TX  C 

C  ;  C 

C . ; . C 

C  ;  C 

C  ;  C 

CALL  P I C 0 U T (  I  BUF2  ,  INAME,  IFLG)  ;  PICBUF  OUTPUT  C 

CALL  RELB(IBUF)  ;  BLOCK  C 

CALL  RELB(IBUF2)  ;  C 

C  ;  C 

C . ; . C 

C  C 

ST0P"<7>FINISHED" 

END 


n  n 


rfc************************************************************** 


c  c 

C  PROGRAM:  LINES. FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN-  15  DEC  8-  SYSTEM:  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  PROGRAM  WILL  FIND  THE  LINES  OF  PRINT  ON  A  C 

C  PAGE  OF  DIGITIZED  TEXT.  C 

C  C 

C  NOTE:  A  BLANK  HORIZONTAL  LINE  MUST  DIVIDE  EACH  C 

C  LINE  OF  PRINTED  TEXT.  C 

C  C 

C  C 

£***«*<.*********************************************************£ 

c  c 

C  C 


PARAMETER  BLEV-16.0D0 
IMPLICIT  INTEGER ( C  ,  D  ,  E  ) 

IMPLICIT  DOUBLE  P8EC I S I  ON ( A ,  B  , T) 

INTEGER  IARRAY1  5  1  2  )  ,  I B  U  F  <  300)  ,  INAME(20) 
INTEGER  D  C  N  T  R (  20)  , ECNTR (20) 

DOUBLE  PRECISION  A  V E (  2  5  8  )  , T A  V E 


c  c 

c . c 

C  ;  C 

ACCEPT "ENTER  FILENAME:"  ;  C 

R  E  A  D (  1  1  ,  1  9  9  )  I N  A  M  E (  1  )  ;  INPUT  FILENAME  C 

19  9  FORMAT ( S  4O )  ;  C 

C.  ;  C 

»  : . ; . C 

G  ;  C 

CALL  P I CFMT (  I  BUF  ,  I  NAME ,  I  FLC  )  ;  C 

CALL  GFMT ( I BUF , NROWS , NCOLS , NB I TS , IMODE )  ;  PICBUF  INPUT  C 

;  BLOCK  C 

CALL  MAKB(IBUF)  ;  C 

CALL  P I C  I  N (  I  BUF  ,  I  NAME  ,  I FLG  )  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

TYPE "LINES  WORKING"  ;  I/O  C 

C  ;  C 

C . ; . C 

G  ;  C 

DO  5  1*1,258  ;  C 

A  V  E  (  I  )  »  0 . 0  D  0  ;  C 

5  CONTINUE  ;  I  N  I T I  A  L I Z  E  C 

;  ALL  C 

DO  1 0  1=1,20  ;  VARIABLES  C 

!)CNTR(  I  )*0  ;  C 

E  C  N  T  R  (  I  )  =  1)  ;  C 

I  II  CONTINUE  ;  C 

C  ;  C 

G . ; . C 
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n  n 


I  • 


c 

c 

c 

c 

c 

c 


•4  0 

3  3 
c 

c .  . 
c 


96 


9  <4 


C 

9  i 
r>  0 


DO  15  J  =  I  ,  I  LN  ;  C 

IF(IARY(J).EQ.l5)GOTO  92  ;  GET  LENGTH  C 

I D= I D+ 1  ;  OF  ROW  C 

CONTINUE  ;  C 

;  C 

I  T=*  IT+  l  ;  C 

IVAL(  IT) *  I D  ;  C 

CONTINUE  ;  C 

;  C 

. ; . C 

;  C 

WRITE  FREE(12)"IT",IT  ;  I/O  C 

;  C 

. ; . C 

C 

JD- IT- 1  ;  C 

DO  35  1=1, JD  ;  C 

J  =*  I  T-  I  ;  C 

DO  40  K= 1 , J  ;  C 

I F (  I VAL ( K )  .  LE .  I  V AL ( K  + 1  )  ) GOTO  90  ;  BUBBLE  SORT  C 

NT’IVAL(K)  ;  ROW  LENGTHS  C 

IVAL(K)«IVAL(K+1 )  ;  C 

I VAL( K+ 1 )=NT  ;  C 

CONTINUE  ;  C 

CONTINUE  ;  C 

;  C 

. ; . C 

;  C 

T  E  M  P  =  0  ;  C 

I  C  =  I  ;  C 

;  C 

N  C  N  T  -  0  ;  C 

DO  43  ICNT*1,IT  ;  C 

I F (  IVAL(ICNT)  .EQ.  I VAL (  IC))NCNT“NCNT+1  ;  C 

CONTINUE  ;  DETERMINE  C 

;  THICKNESS  OF  C 

I  F  (  NCNT  .  LE  .  TEMP  )  GOTO  9-,  ;  ROW  LENGTHS  C 

TEMP-NCNT  ;  C 

I T  E  M  P  =  1 0  ;  C 

;  C 

;  C 

I C- I C+ 1  ;  C 

IF(IC.GT.IT)GOTO  95  ;  C 

GOTO  96  ;  C 

. ; . C 

;  C 

DO  )0  IH=1,IT  ;  C 

WRITE  FREE(L2)”IVAL",IH,IVAL(IH)  ;  C 

CONTINUE  ;  C 

;  I/O  C 

MAX ■ I  V AL (  I TEMP )  ;  C 

;  C 

WRITE  FREE( 1 2 ) "MAX" , MAX  ;  C 

;  C 

. ; . C 


I  I  2 


c 

PROGRAM:  RL.MX.FR  LANGUAGE:  FORTRAN  V  C 

C 

i:  WRITTEN:  I  )  DEC  8  •*  SYSTEM:  ECLIPSE  C 

C 

C  PURPOSE  :  THIS  PROGRAM  DETERMINES  THE  THICKNESS  OF  A  LETTER.  C 
G  C 

f***************************************************************^ 

C  C 

INTEGER  I  BUF (  300  ) 

INTEGER  I  ARY (  2  5  6  ) 

INTEGER  INAME ( 20 ) 

INTEGER  IVAL(IOO) 

C . C 

C  ;  C 

ACCEPT "ENTER  FILENAME:"  ;  C 

READ(  l  I  ,  199)  I  NAME (  1  )  ;  ENTER  FILENAME  C 

199  FORMAT ( S40 )  ;  C 

C . ; . C 

C  ;  C 

ACCEPT" ENTER  TOP  ROW:  ".1ST  ;  ENTER  C 

ACCEPT" ENTER  BOTTOM  ROW;  ",LST  ;  WINDOW  C 

ACCEPT" ENTER  LEFT  COLUMN:  ",ICO  ;  COORDINATES  C 

ACCEPT" ENTER  RIGHT  COLUMN:  ",IRC  ;  C 

C . ; . C 

C  ;  C 

CALL  P I CFMT(  I  BUF  ,  I  NAME ,  I FLG )  ;  PICBUF  INPUT  C 

CALL  GFMT (  I BUF , NROWS , NCOLS  , NB ITS ,  IMODE  )  ;  BLOCK  C 

CALL  MAKB(IBUF)  ;  C 

CALL  PICI.N(  I  BUF  ,  INAME  ,  IFLG)  ;  C 

C . ; . C 

C  ;  C 

C  ;  C 

I LN= I RC- I C0  +  1  ;  C 

I T  *  0  ;  INITIALIZE  C 

;  ALL  C 

DO  1  JNM.100  ;  VARIABLES  C 

I  V  A  L  (  J  N  >*99999  ;  C 

1  CONTINUE  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

C  ;  C 

DO  5  I R= I  ST  ,  LST  ;  C 

CALL  CROW(  I BUF  ,  I R  ,  ICO  ,  I LN  ,  I  AR Y  )  ;  C 

;  C 

DO  10  I  *  1  ,  I  LN  ;  C 

IF( IARY( I ) . EQ. 1 5 )GOTO  10  ;  GET  TO  FIRST  C 

GOTO  91  ;  BLACK  PIXEL  C 

10  CONTINUE  ;  ONLINE  C 

;  C 

GOTO  5  ;  C 

;  C 

9  1  I  D  -  0  ;  C 


1  1  1 


CALL  PICOUT(  IBUF2  ,ONAME,  IFLG) 
CAL.L  RELB (  I BUF  ) 

CALL  RELB(  IBUF2  ) 

STOP" <7  >F I N I  SHED" 

END 


1  1  0 


DO  5  I -YR  <  l  )  , YR( 3  ) 

CALL  GROW(  I BUF  ,  I  ,  1  , NCOLS  ,  I  ARRAY ) 
DO  10  J  *  l  ,  2  5  6 

I  ARRAY ( J ) *NOT (  IARRAY (  J  )  ) 

10  CONTINUE 

CALL  PROW(  I BUF  ,  I  ,  1  , NCOLS  ,  I  ARRAY  ) 
5  CONTINUE 

TY PE " N EGAT I VE  FORMED" 


C  THIS  PROGRAM  WILL  DIVIDE  A  PICTURE  INTO  FOUR  QUADRANTS 

C  AND  INTERCHANGE  EACH  OF  ITS  DIAGONAL  QUADRANTS.  YOU 

C  MUST  PROVIDE  THE  UPPER  LEFT  HAND  AND  THE  LOWER  RIGHT 

C  HAND  COORDINATES  OF  THE  ORIGINAL  PICTURE. 

XCT-  ( XC (  3 ) +XC (  l  )  )  /  2 
XCT- I  NT ( XCT ) 

XC ( 2 ) -XCT  + 1 


YR ( 2 ) J ( YR(  3  )  +  Y  R  (  1  )  )/2 
YR( 2 )  =  INT( YR( 2  )  ) 


I C  L  N (  1  )«XC( 2 )-XC<  1  ) 
ICLN ( 2 ) -XC ( 3 ) -XCT 


IC(  1  )-  1 

IC(  2  ) -2  5  7  - ICLN(  1  ) 


I  RW(  1  ) - YR ( 2 ) - YR (  1  )  +  1 
IRC  1  )  -  2  5  7  -  I  R  W  (  1  ) 

I  R  (  2  )  -  1 


DO  15  J  -  l  ,  2 

KR- I R ( J ) 

L-  1 

DO  20  IROW-YRC J ) , YR ( J+l ) 

CALL  GROW(  IBUF  ,  IROW,XC(L)  ,  ICLN(L)  ,  IARRAY) 

CALL  PROW(  IBUF2  ,KR,  IC ( L+ 1  )  ,  ICLN(L)  ,  IARRAY) 

CALL  GROW(  IBUF,  I  ROW , XC( L+ 1  )  ,  ICLNCL+l  )  ,  IARRAY) 
CALL  PROW(  IBUF2,KR,IC(L)  , ICLN ( L+ 1  )  .IARRAY) 

KR-K.R  +  l 

20  CONTINUE 

Y  R ( 2 )-YR( 2 )  +  l 
l 5  CONT I NUE 
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£***************************************************************£ 

C  C 

c  c 

C  ORIGINAL  PROGRAM:  NEG.FR  BY  R.  A.  SIMMONS  C 

c  c 

c  c 

C  PROGRAM:  QSHIFT.FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  PROGRAM  WILL  PERFORM  THE  FOLLOWING  OPERATIONS :C 
C  1.  CREATE  THE  NEGATIVE  OF  A  SECTION  OF  A  VIDEO  FILE  C 

C  2.  DIVIDE  THE  SECTION  INTO  FOUR  QUADRANTS  C 

C  3.  INTERCHANGE  EACH  DIAGONAL  QUADRANT  C 

C  C 

£***************************************************************£ 
C  C 

c  c 


IMPLICIT  I NTEG  E  R ( X , Y  ,  Z ) 

INTEGER  I  A  R  R  A  Y (  300)  ,  IBUF(300)  ,  IBUF2(300) 

INTEGER  INAME ( 20 ) , ONAME ( 20 ) 

INTEGER  XC(4)  ,YR(4)  ,  ICLN( 2  ) 

INTEGER  IC( 2 )  ,  IR( 2 )  ,  IRW< 2) 

C  C 

C  C 

ACCEPT "ENTER  INPUT  FILENAME:  " 

R  E A  D (  11,199)  I N  AME (  1  ) 

ACCEPT" ENTER  OUTPUT  FILENAME:  " 

R  E  A  D (  1  1  ,  19  9) ONAME (  1  ) 

199  FORMAT  ( S  4  0  ) 


991  TYPE "ENTER  ROW  AND  COLUMN  OF  CORNERS:" 

ACCEPT” UPPER  LEFT  HAND  CORNER  - - > "  ,  Y R (  1  )  , X C (  1  ) 
ACCEPT" LOWER  RIGHT  HAND  CO R N E R - - > "  ,  Y R ( 3  )  ,  X C ( 3  ) 


I F ( (XC( 1 ) . GE . XC ( 3) ) .AND. (YR( 1 ) .GE. YR( 3) ) )GOTO  991 

C  ORIGINAL  PROGRAM:  NEG.FR  BY  R.  A.  SIMMONS 
C  REWRITTEN:  NEG1.FR 

C  THIS  PROGRAM  WILL  CREATE  THE  NEGATIVE  OF  A  VIDEO  FILE. 

TYPE"QSHIFT  WORKING" 

CALL  PICFMT(  IBUF  ,  INAME  ,  IFLG) 

CALL  GFMT(IBUF,NROWS,NCOLS,NBITS,  I  MODE) 

CALL  M A  K  B (  IBUF) 

CALL  P IC I N (  I BUF ,  I  NAME ,  I FLG ) 

CALL  GFMT (  IBUF, NROWS.N COLS, NBITS,  I  MODE) 

CALL  P  FMT (  IBUF2,NROWS,NCOLS, NBITS,  IMODE) 

CALL  MAKB (  I  BUF2  ) 
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DO  20  IR-IST.LST 

CALL  GROW(  I  BUF  ,  IR ,  1  ,  2  5  6  ,  IARY ) 

DO  25  J  «  l  ,256 

I F ( I AR Y ( J ) . EQ . 1 5 )GOTO  25 
CALL  P  P  NT  (  IBUF2,  IR,  J,  1  ) 
CONTINUE 

CONTINUE 


CONVERT 


(0/15) 

GREYSCALE 

TO 

(1/0) 

GREYSCALE 


DO  3  5  I R- I  ST , LST 

CALL  GROW(  IBUF2  ,  IR ,  1  ,  256  ,  IARY) 

DO  40  J-2  ,  2  5  5 

I  F (  I  ARY ( J  )  .  EQ . 0 )G0T0  40 

CALL  GROW( I BUF2 , IR- 1 , J- 1 , 3 , ICRY ) 
CALL  GROW(  I BUF2  ,  I  R  ,  J-  1  ,  3  ,  I  DR Y ) 
CALL  GROW( I BUF2 , I R+l , J- 1 , 3 , I ERY ) 
IX-0 

DO  4  5  K *  1  ,  3 

IX»ICRY(K)+IDRY(K)+IERY(K)+IX 

CONTINUE 

CALL  PPNT(  I  BUF3  ,  IR , J  ,  IX ) 

CONTINUE 

CONTINUE 


CALCULATE 
NUMBER  OF 
NEAREST 
NEIGHBORS 


TYPE" < 7 >MTRXA  FINISHED" 
RETURN 
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PROGRAM:  MTRXA.FR 


LANGUAGE:  FORTRAN  V 


WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE 

PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  SPOTS. FR  AND 

VTHIN.FR  .  IT  CALCULATES  MATRIX  A  WHERE  EACH 
ELEMENT  IS  THE  NUMBER  OF  NEAREST  NEIGHBORS 
THAT  A  PIXEL  HAS  (MAXIMUM  IS  9,  MINIMUM  IS  0) 


******************************************************* 


SUBROUTINE  MTRXA 


COMMON  I BUF (  300  ) 
COMMON  I BUF2 ( 300) 
COMMON  I  B  U  F  3  ( 300) 
COMMON  IARY(256) 
COMMON  I  NAME (  20  ) 
COMMON  IFLAG 
COMMON  1ST 
COMMON  LST 
RE^AL  WFILE(  5) 
INTEGER  I C  R  Y ( 4  ) 
INTEGER  I D  R  Y (  4  ) 
INTEGER  IERYC4) 


TYPE" MT  R  X  A  WORKING 


DO  5  I  »  l  ,  2  5  6 

I AR Y (  I  )  -0 
CONTINUE 

I  X-0 


I/O 


INITIALIZE 

ALL 

VARIABLES 


WFILE="BK" 


CALL  P I C  I  N (  I  BUF2  , WF I LE , 0  ) 
CALL  P I C I N ( I BUF3 , WF I LE , 0 ) 


TYPE "CHANGED  TARGET" 


CALL  PICOUT(  IBUF2  ,  INM2  ,  IFLG) 


GOTO  91 


ER»E 1 / E2 

WRITE  FREE (  1  0 ) "TEMPLATE  IS",ER, 
/  "  TIMES  THE  TARGET  ENERGY" 


I  F ( ER . GE .  1  .  5  DO  )GOTO  77 
TYPE "ENERGIES  ARE  THE  SAME" 
GOTO  91 


ENERG I ES 
ARE 

THE  SAME 


ER-E2 / E 1 

E3»IDINT((ER*15.0D0)+0.5D0) 

E4-E3+1 

IE3D-IABS((IE1*15)-(IE2*E3)) 

IE4D=IABS((IE1*15)-(IE2*E4)) 

I  F (  IE3D.GE.  IE4D)E3*E4 

IF(E3.LT.O)E3-0 
I F ( E  3 . GT  .  1 5 ) E3- 1 5 

DO  57  I  R  *  3  9 , 6  5 

CALL  GROW(  I  BUF 1  ,  I R ,  I  CL ,  I LN ,  I  ARY ) 

DO  38  J»1 , I LN 

I F ( IARY(J) .EQ.0)IARY(J)“1 5-E3 
CONTINUE 

CALL  PROW(IBUFl  ,  I R ,  I  CL ,  I LN ,  I  ARY ) 
CONTINUE 

TYPE "TEMPLATE  CHANGED" 

CALL  P I COUT(  I  BUF  1  ,  I NM 1  ,  I FLG ) 


TEMPLATE  ENERGY 
IS  GREATER  THAN 
THE  WINDOW 
ENERCY 


SET  TEMPLATE 
ENERGY  EQUAL 
TO  WINDOW 
ENERGY 

SET  THE 
TEMPLATE 
ENERGY  EQUAL 
TO  THE  WINDOW 
ENERGY 
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O  O 


c . c 

C  ;  c 

C  ;  C 

WRITE  FREE (  1 0 ) "TEMPLATE  ENERGY", El  ;  I/O  C 

WRITE  FREE (  1 0 ) "TARGET  ENERGY"  ,  E2  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

C  ;  DETERMINE  C 

C  ;  IF  WINDOW  ENERGY  C 

I F ( E2  .  LT . E 1  )GOTO  80  ;  IS  GREATER  THAN  C 

C  ;  THE  TEMPLATE  ENERGY  C 

C  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

C  ;  C 

ER*E2/E1  ;  C 

WRITE  FREE ( 1 0 ) "TARGET  IS",ER,  ;  I/O  C 

/  "  TIMES  THE  TEMPLATE  ENERGY"  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

C  ;  DETERMINE  IF  C 

I F( ER . GE .  1  . 0  1  DO  )COTO  67  ;  THE  ENERGIES  C 

TYPE "ENERGIES  ARE  THE  SAME"  ;  ARE  EQUAL  C 

GOTO  91  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

C  ;  C 

67  ER-E1/E2  ;  C 

E3* ID INT( ( ER* 1 5 . ODO ) +0 . 5D0  )  ;  C 

E4-E3+1  ;  C 

;  C 

;  WINDOW  ENERGY  C 

I  E3D- I ABS <  (  IE l M 5 ) - <  IE2*E3 )  >  ;  IS  GREATER  C 

I E4D- I ABS (  ( IE  1  *  1  5 ) - (  I E2 *E4 )  )  ;  THAN  THE  C 

I F ( IE3D.GE. IE4D)E3-E4  ;  TEMPLATE  ENERGY  C 

;  C 

;  C 

I F ( E3 . LT . 0 ) E3-0  ;  C 

IF<  E3 .GT.  1 5  )  E  3  *  15  ;  C 

;  C 

;  C 

DO  40  IR-IST.LST  ;  SET  THE  C 

CALL  GROW( IBUF2 , IR , ICO, ILX, IARY)  ;  WINDOW  ENERGY  C 

;  EQUAL  TO  C 

DO  45  J-l.ILN  ;  THE  TEMPLATE  C 

I F (  I AR Y ( J )  . EQ . 0 ) I  ARY ( J ) -  1 5-E3  ;  ENERGY  C 

*5  CONTINUE  ;  C 

;  C 

CALL  PROW(  I BUF2  ,  I R ,  ICO ,  I  LX  ,  I  ARY  )  ;  C 

;  c 

hO  CONTINUE  ;  C 
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C . ; . C 

C  ;  C 

E l -0 . 000  ;  C 

E2-0.0D0  ;  INITIALIZE  C 

E3-0  ;  ALL  VARIABLES  C 

E  4  =  0  ;  C 

C . ; . C 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


C  ;  C 

ICL-ICOL(INUM)  ;  C 

ILN=ILEN< INUM)  ;  C 

DO  5  IR-39,65  ;  C 

CALL  GROW(  I BUF 1  ,  I R ,  ICL ,  ILN ,  I  ARY )  ;  C 

;  CALCULATE  C 

DO  l 0  J- 1 , I LN  ;  ENERGY  C 

IF( IARY( J) .EQ.0)E1-E1+1 .ODO  ;  IN  TEMPLATE  C 

10  CONTINUE  ;  C 

;  c 

5  CONTINUE  ;  C 

;  C 

E 1 »DSQRT( E I )  ;  C 

I E 1  *  I DINT( E 1 )  ;  C 

C . ; . C 

C  ;  C 

I  LX- I LC- ICO+ L  ;  C 

DO  25  IR-IST.LST  ;  C 

CALL  GROW(  I BUF2  ,  I R ,  ICO ,  I  LX ,  I  AR Y  )  ;  C 

;  C 

DO  30  J- L , ILX  ;  C 

I F ( I AR Y ( J ) . EQ . 0 ) E2 =E2+ 1 . ODO  ;  CALCULATE  C 

30  CONTINUE  ;  ENERGY  C 

;  IN  WINDOW  C 

25  CONTINUE  ;  C 

;  C 

E2-DSQRT(E2)  ;  C 

IE2-IDINT(E2  )  ;  C 

C . ; . C 


1  02 


£**********************★********★**★****************************£ 

C  C 

c  c 

C  PROGRAM:  MEND.FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  PROGRAM  WILL  MODIFY  THE  ENERGY  OF  A  LOW  C 

C  ENERGY  INPUT  TO  EQUAL  A  HIGH  ENERGY  INPUT.  C 

C  NOTE:  THIS  PROGRAM  READS  FROM  A  STANDARD  TEMPLATE  FILE  C 

C  (FIGURE  7.1). FOR  RESEARCH,  THE  WINDOW  ENERGY  WAS  C 

C  ALWAYS  SET  TO  THE  TEMPLATE  ENERGY.  TEMPLATE  C 

C  CHARACTERS  ALL  HAD  ABOUT  THE  SAME  ENERGY.  C 

C  C 

£***************************************************************£ 
C  C 


INTEGER  I NM 1 ( 7 ) , I NM2 ( 7 ) 

INTEGER  IARY(256)  ,  ICOL(8)  ,  ILEN(8) 
INTEGER  I BUF l (  300  )  , I BUF2 (  300 ) 
INTEGER  I NUM , E  3 , E  4 
DOUBLE  PRECISION  E N 1  , E N 2  , E R  ,  E D 


C . C 

C  ;  C 

ACCEPT" ENTER  TEMPLATE  FILENAME:  "  ;  C 

READ ( l 1 , 1 99 ) INM 1 ( 1 )  ;  ENTER  C 

;  FILE  NAME  C 

ACCEPT" ENTER  TARGET  FILENAME:  "  ;  C 

R  E  A  D (  11,199)  I N  M  2 (  1  )  ;  C 

199  FORMAT ( S  20 )  ;  C 

C . ; . C 

C  ;  C 

TYPE "ENTER  TEMPLATE  TO  MEND:  "  ;  ENTER  C 

TYPE"1->A,  2  - >  B ,  3  - >  C ,  4 - >  D "  ;  TEMPLATE  C 

TYPE"5->E,  6->F,  7->G,  8->H"  ;  C 

ACCEPT" ENTER  NUMBER:  ",INUM  ;  C 

C  .  . . ; . C 

C  ;  C 

ACCEPT" ENTER  TOP  ROW:  ".1ST  ;  C 

ACCEPT" ENTER  BOTTOM  ROW:  ",LST  ;  ENTER  WINDOW  C 

ACCEPT "ENTER  LEFT  COLUMN:  ",ICO  ;  COORDINATES  C 

ACCEPT" ENTER  RIGHT  COLUMN:  ",ILC  ;  C 

C . ; . C 


C 

c 


INM1 -"TEMPSET . OR"  ;  C 

CALL  P ICFMT (  I  BUF  1  ,  INM 1  ,  I FLG  )  ;  C 

CALL  PICFMTC  IBUF2  ,  INM2  ,  IFLG)  ;  C 

;  c 

CALL  GFMT(  I BUF 1  , NROWS , NCOLS  ,  NB  ITS  ,  IMODE )  ;  PICBUF  INPUT  C 

;  BLOCK  C 

CALL  M A  K  B (  I  BUF  1  )  ;  C 

CALL  MAK  B (  I BUF2  )  ;  C 

;  C 

CALL  P I C I N (  I  BUF  1  ,  I NM 1  ,  I FLG  )  ;  C 

CALL  P I C I N (  I BUF2  ,  I NM2  ,  I FLG  )  ;  C 


1  0  1 


DO  7  5  I  =  1  ,  C-  1 

TYPE "LINE"  ,  I  ,  "  IS  AT  ROW"  , DCNTR (  I  ) 
ECNTR (  I  ) 

CONTINUE 


OPEN  I , "LDAT.TX" 

WRITE  B I  NARY (  l  ) DCNTR  ,  ECNTR  ,  C 

TYPE"<7>FINISHED  " 


TO"  , 


C 

C 


I/O 


STORE  LINE 
POSITIONS  IN 
LDAT . TX 


I/O 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


CALL  RELB(IBUF) 


PICBUF 

OUTPUT 


O  O 


c . 

c 

DO  20  1*1,240  ;  C 

CALL  GROW(  I BUF  ,  I  ,  l  , NCOLS ,  I  ARRAY )  ;  C 

DO  25  J  -  1  ,  2  56  ;  CALCULATE  C 

K- I  ARRAY ( J )  +  1  ;  AVERAGE  PIXEL  C 

AVE (  1  +  1  ) - AVE (  1  +  1  )  +  DFLOAT( K )  ;  VALUE  FOR  C 

25  CONTINUE  ;  EACH  ROW  C 

AVE (  1  +  I  ) -AVE(  1  +  I  ) / 2 5 6 . 0D0  ;  C 

20  CONTINUE  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

AVE ( 1 ) = 1 6 . ODO  ;  C 

DO  30  1-242,258  ;  SET  TOP  C 

AVE(  I ) -  1 6  .  ODO  ;  AND  BOTTOM  TO  C 

30  CONTINUE  ;  WHITE  C 

C . ;  .  .  .  . . C 

C  ;  C 

TAVE-0.0D0  ;  CALCULATE  C 

DO  40  1-2,241  ;  AVERAGE  PIXEL  C 

TAVE-TAVE  +  AVE (  I  )  ;  VALUE  FOR  C 

40  CONTINUE  ;  ENTIRE  VIDEO  C 

TAVE-TAVE/256 . ODO  ;  FILE  C 

C . ; . C 

C  ;  C 

WRITE  FREE (  1 0 ) "TAVE-  "  , T  A  V  E  ;  C 

I F( TAVE . GT . 1 0 )GOTO  991  ;  DETERMINE  IF  C 

TYPE "FILE  CANNOT  BE  PROCESSED"  ;  FILE  CAN  C 

GOTO  9999  ;  BE  PROCESSED  C 

991  TYPE "FILE  MAY  BE  PROCESSED"  ;  C 

C . ; . C 

C  ;  C 

C-  l  ;  C 

DO  451-2,  242  ;  C 

I F (  ( AVE (  I  )  . NE . BLEV )  ;  DETERMINE  C 

/  .AND.  ;  LOWER  EDGE  OF  C 

/  ( AVE (  I -  1  )  . EQ .  BLEV )  )GOTO  91  ;  LINE  C 

GOTO  45  ;  C 

9  1  DCNTR ( C )  -  I  -  l  ;  C 

C-C+l  ;  C 

45  CONTINUE  ;  C 

C . ; . C 

C  ;  C 

C-  1  ;  C 

DO  50  1-2  ,  25  7  ;  C 

I F (  ( AVE (  I  )  .  NE  .  BLEV  )  ;  DETERMINE  C 

/  .AND.  ;  UPPER  EDGE  OF  C 

/  ( A  V  E ( I  +  1  )  . EQ. BLEV) )GOTO  92  ;  LINE  C 

GOTO  50  ;  C 

92  ECNTR ( C )  -  I  ;  C 

C-C+l  ;  C 

50  CONTINUE  ;  C 

C  ;  C 

C . ; . C 
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n  o 


PROGRAM:  SPOTS. FR 


LANGUAGE:  FORTRAN  V 


WRITTEN:  15  DEC  84 


SYSTEM:  ECLIPSE 


PURPOSE:  THIS  PROGRAM  WILL  REMOVE  THE  SPOTS  IN  THE  INPUT.  C 
IT  DELETES  POINTS  OF  ONE  PIXEL,  FILLS  HOLES  OF  C 

ONE  PIXEL,  AND  DELETES  TIPS  OF  SPIKES.  C 

C 
C 

***************************************************£ 

C 

COMMON  IBUF(300) 

COMMON  IBUF2(300) 

COMMON  IBUF3(300) 

COMMON  IARY(256) 

COMMON  INAME(20) 

COMMON  I  FLAG 
COMMON  1ST 
COMMON  LST 
INTEGER  IBRY(256) 

REAL  WF I LE ( 5  ) 


ACCEPT "ENTER  INPUT  FILENAME: 
R  E  A  D (  11,199)  I  NAME (  1  ) 

F ORMAT ( S  4  0 ) 


INPUT  FILENAME 


TYPE" SPOTS  WORKING" 


CALL  PICFMT( IBUF , INAME , IFLG) 

CALL  G  FMT  (  I BUF , NROWS , NCOLS , NB ITS ,  I  MODE) 
CALL  MAKB(IBUF) 

CALL  PICIN(  IBUF,  INAME,  IFLG) 

WF I LE- " BK" 

CALL  G  FMT (  IBUF, NROWS, NCOLS, NBITS,  I  MODE) 
CALL  PFMT(  IBUF2  .NROWS  .NCOLS  , NBITS  ,  IMODE) 
CALL  MAK  B (  I BUF2  ) 

CALL  P I C I N (  IBUF2  ,WFILE,  IFLG) 

CALL  G  FMT (  IBUF, NROWS, NCOLS, NBITS,  IMODE) 
CALL  PFMT(IBUF3, NROWS, NCOLS, NBITS,  I  MODE) 
CALL  MAKB (  I BUF  3 ) 

CALL  P I C I N (  I BUF 3 , WF I LE  ,  I FLG  ) 


PICBUF  INPUT 
BLOCK 

(BK  IS  A  BLACK 
VIDEO  FILE) 


C . ; . C 

C  ;  C 

I  ST*  2  ;  SET  TOP  ROW  C 

L  S  T  *  2  4 1  ;  SET  END  ROW  C 

C . ; . C 

C  ;  C 

CALL  MTRXA  ;  FORM  MATRIX  A  C 

C . ; . C 

C  ;  C 

WRITE  FREE(  1 0 ) " BEG  I N  1"  ;  C 

DO  5  IR-2 , 24 1  ;  C 

CALL  GROW(  I  B  U  F  3  ,  IR,  l  ,256,  IBRY)  ;  C 

DO  1 0  J-2  ,  2  5  5  ;  C 

I F(  I  BRY < J )  . NE . 0 ) GOTO  10  ;  C 

CALL  GPNT( IBUF3  ,  IR-1  ,  J  ,  IARY(  1  )  )  ;  C 

CALL  GPNT (  I  BUF3  ,  I  R  ,  J  -  1  ,  I  ARY ( 2 )  )  ;  C 

CALL  GPNT (  I BUF3 , IR , J+ 1  ,  I  ARY ( 3 )  )  ;  DETERMINE  IF  C 

CALL  GPNT(  I BUF3  ,  IR+ 1  ,  J  ,  I ARY( 4  )  )  ;  POINT  IS  A  C 

DO  1 5  K- 1  ,  4  ;  HOLE  C 

I F (  I  ARY ( K  )  .  LT  .  1  )GOTO  10  ;  C 

15  CONTINUE  ;  C 

CALL  PPNT(  I BUF  ,  I  R , J , 0 )  ;  C 

10  CONTINUE  ;  C 

5  CONTINUE  ;  C 

WRITE  FREE( 10) "END  1"  ;  C 

C . . ; . C 

C  ;  C 

WRITE  FREE (  l 0 ) " BEG  I N  2"  ;  C 

DO  30  IR-2, 241  ;  C 

CALL  GROW( I BUF3  ,  IR  ,  1  ,  2  5  6  ,  I  BRY  )  ;  C 

DO  35  J-2,255  ;  DETERMINE  IF  C 

IF(  IBRY(J)  .NE.  1  )GOTO  35  ;  POINT  IS  A  C 

CALL  PPNT( I BUF  ,  IR  ,  J ,  1 5 )  ;  SPOT  C 

35  CONTINUE  ;  C 

30  CONTINUE  ;  C 

WRITE  FREE(lO)"END  2"  j  C 

C . . . ; . C 


C  ;  C 

WRITE  FREE( 10) "BEGIN  3"  ;  C 

DO  50  IR-2,  241  ;  C 

CALL  GROW(  I BUF3  ,  IR ,  1  ,  2  5  6  ,  I BRY  )  ;  C 

DO  55  J-2,255  ;  C 

I F ( I BRY ( J  )  . NE .  2 )GOTO  55  ;  C 

CALL  GPNT(  IBUF3  ,  IR-l  , J-l  ,  IARY(  1  )  )  ;  C 

CALL  GPNT(  I BUF3  ,  IR- 1  ,  J , I  ARY ( 2  )  )  ;  C 

CALL  GPNT(  IBUF3  ,  IR- 1  , J+l  ,  IARY( 3)  )  ;  C 

CALL  GPNT(  I BUF3  ,  IR  ,  J- 1  ,  I  ARY ( 4  )  )  ;  DETERMINE  C 

CALL  GPNT<  I  BUF3  ,  I R , J+ 1  ,  I  ARY (  5 )  )  ;  IF  C 

CALL  GPNT (  I  BUF3  ,  I R+ 1  , J- 1  ,  I  ARY ( 6  )  )  ;  POINTS  C 

CALL  GPNT(  IBUF3  ,  IR+l  , J  ,  IARY( 7)  )  ;  ARE  SPOTS  C 

CALL  GPNT( I BUF3  ,  IR+l  , J+l  , IARY( 8  )  )  ;  AND  TIPS  C 

;  OF  SPIKES  C 
;  C 
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(«r 


60 


65 

9  1 

5  5 
5  0 


C 

C 


90 

94 

85 

80 

C  .  .  .  . 

C 

9999 


C 

C 


DO  60  K= 1  ,  8  ;  C 

IF(  I AR  Y ( K )  .  EQ. 2 )GOTO  91  ;  C 

CONTINUE  ;  C 

;  C 

DO  65  L»l,8  ;  C 

IF(  IARV(L)  .  EQ.  5 )GOTO  91  ;  C 

I F( I AR Y ( L ) . EQ . 6 )GOTO  91  ;  C 

IF( IARY( L)  . EQ. 7 )GOTO  91  ;  C 

CONTINUE  ;  C 

GOTO  55  ;  C 

CALL  PPNT(  I  BUF  ,  I  R  ,  J  ,  1  5  )  ;  C 

CONTINUE  j  C 

CONTINUE  ;  C 

;  C 

WRITE  FREE( 10) "END  3"  ;  C 

. ; . C 

C 

WRITE  FREE( 1 0 ) "BEGIN  4"  ;  C 

;  C 

DO  80  IR-2  ,  24  1  ;  C 

CALL  GROW( I BUF3 , IR , 1 , 2 5 6 , I BRY )  ;  C 

DO  8  5  J»2 , 2  5  5  ;  C 

I  F (  I BRY ( J )  . NE . 4 )G0T0  8  5  ;  C 

CALL  GPNT(  IBUF3 ,  IR- 1  , J- 1  ,  IARY(  1  )  )  ;  C 

CALL  GPNT( I BUF3 , IR- 1  , J  ,  J ,  I  ARY ( 2 )  )  ;  C 

CALL  GPNT( I BUF3 , IR- 1 , J+l , I ARY( 3 ) )  ;  C 

CALL  GPNT( I BUF3 , I R , J- 1  ,  I  ARY ( 4 ) )  ;  C 

CALL  GTNT( I BUF3 , IR , J  + 1 , I ARY( 5 ) )  ;  DETERMINE  IF  C 

CALL  GPNT( I BUF3  ,  IR+1  ,  J-  1  ,  I ARY( 6  )  )  ;  POINTS  ARE  C 

CALL  GPNT(  IBUF3  ,  IR+1  , J,  IARY(7)  )  ;  TIPS  OF  C 

CALL  GPNT(  IBUF3,  IR+1  , J+l  ,  IARY(8)  )  ;  SPIKES  C 

;  C 

I BRY (  1  ) - I ARY(  1  )  +  I ARY( 2 )  +  I ARY ( 3  )  ;  C 

I  BRY  (  2  )  <■  I  ARY  (  1  )  +  I  ARY  (  4  )  +  I  ARY  (  6  )  ;  C 

I BRY ( 3 ) - I  ARY ( 3 )  +  I ARY ( 5 )  +  I ARY ( 8  )  ;  C 

I BRY ( 4 ) - I  ARY ( 6 )  +  I AR Y ( 7 )  +  I ARY ( 8  )  ;  C 

;  C 

DO  90  K»l,8  ;  C 

I F( I BRY ( K) . EQ . 0 )GOTO  94  ;  C 

CONTINUE  ;  C 

GOTO  85  ;  C 

CALL  PPNT(  I  BUF ,  IR  ,  J ,  1 5 )  ;  C 

CONTINUE  ;  C 

CONTINUE  ;  C 

. ; . C 

;  c 

CALL  PICOUT( IBUF , INAME , IFLG)  ;  PICBUF  C 

CALL  RELB(IBUF)  OUTPUT  C 

CALL  RELB(IBUF2)  ;  BLOCK  C 

CALL  R  E  L  B (  I BUF3  )  ;  C 

. ; . C 

C 

ST0P"<7>SP0TS  FINISHED" 

END 
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* 


CALL  GPNT( I BUF3 ,  I I - 1  ,  J  J  -  1  ,  I VAL(  1  )  ) 
CALL  GPNT(IBUF3,II-1,JJ+1,IVAL(2)) 
CALL  GPNT( I  B  U  F  3 , I  I  +  1 ,JJ-1 , I VAL( 3 ) ) 
CALL  GPNT(  IBUF3  ,  I  I  +  1  , J J+ 1  ,  IVAL( 4  )  ) 

DO  7  1-1,4 

IF( I VAL( I ) . EQ . 8 >GOTO  100 
7  CONTINUE 


CALL  GPNT(  I8UF3  ,  I  I  -  1  , J J , I VAL( 5  )  ) 

CALL  GPNT (  IBUF3  ,  II  ,  J  J -  1  ,  IVAL(6)  ) 

CALL  GPNT(  IBUF3  ,  I  I  ,  JJ+1  ,  I VAL( 7 )  ) 

CALL  GPNT(  I BUF3  ,  I  1  +  1  , J J ,  I VAL( 8 )  ) 

I F ( (  I VAL(  1  )  . NE . 6 )  .OR.  ( I  V A L (  1  )  .  N E  .  7  )  )  G OTO  91 
IF(  (  I  V  A  L ( 5 )  .NE.O)  .OR.  (  I  V A L ( 6 )  .  N E  .  0  )  )  G OTO  100 


91  I F (  ( I VAL( 2 )  . NE . 6 )  .OR.  ( I  V A L ( 2  )  . N E . 7 ) ) GOTO  92 

IF( (  IVAL( 5)  .NE.O)  .OR.  (  I  V A L ( 7  )  . N E  .  0  )  ) GOTO  100 

92  IF(  (  IVAL( 3)  . NE.  6)  .OR.  ( I  V A L ( 3 )  .  N E  .  7  )  )  G OTO  93 

IF<  (  IVAL(6)  .NE.O)  .OR.  (  I  V A L ( 8  )  .  N E  .  0  )  ) G OTO  100 

93  I F (  (  I VAL ( 4 )  . NE . 6 )  .OR.  ( I  V A L ( 4  )  . N E  .  7  )  ) GOTO  94 

IF(  (  IVAL(  7)  .NE.O)  .OR.  (  I  V A L ( 8  )  .  N E  .  0  )  ) GOTO  100 


GOTO  94 

100  CALL  PPNT(  I  BUF  ,  I  I  , JJ  ,  1  5  ) 

94  RETURN 

END 
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PROGRAM:  TDT4.FR 


LANGUAGE-  FORTRAN  V 


WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE 


PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  THIN.FR  TO  TEST  FOR 
D  2  AND  NEIGHBOR  VALUE  OF  9  DELETION. 

********************************************************* 

SUBROUTINE  TDT4(II,JJ) 

COMMON  IBUF( 300) 

COMMON  IBUF2(300) 

COMMON  IBUF3(300) 

COMMON  IARY(256) 

COMMON  I  NAME (20) 

COMMON  IFLAG 
COMMON  1ST 
COMMON  LST 
INTEGER  I  V  A  L ( 9  ) 

DO  5  1-1,9 

I  V  A  L (  I  )  -0 
CONTINUE 

CALL  GPNT(  IBUF3,  II-l  , J  J -  1  ,  I  V  A  L (  1  )  ) 

CALL  G  P  N  T (  IBUF3,  I  I -  1  ,  J  J  +  1  ,  IVAL(2)) 

CALL  GPNT(  IBUF3 ,  II  +  l  ,  J  J -  1  ,  I  V  A  L ( 3 )  ) 

CALL  GPNT(  IBUF3  ,  II  +  l  ,  J  J  +  l  ,  I VA  L( 4  )  ) 

DO  7  1-1,4 

I F ( I VAL( I ) . EQ . 9 )GOTO  100 
CONTINUE 

CALL  GPNT(  I BUF3  ,  I  I  -  1  , J J ,  I VAL(  5  )  ) 

CALL  GPNT (  IBUF3,  I  I  ,  J  J  -  l  ,  I  VAL ( 6  )  ) 

CALL  GPNT(  IBUF3  ,  II  ,  J  J  +  1  ,  I  V  A  L ( 7 )  ) 

CALL  GPNT( IBUF3, II+l ,JJ,IVAL(8)) 

I F (  (  I VAL(  1  )  .  LE .  3 )  .OR.  (  I  V A L (  1  )  . G E  .  9  )  )  GOTO  91 
I F (  (  I VAL(  1  )  . EQ . 4  ) 

/  . OR .  (  I  VAL(  1  )  .  EQ  .  6  ) 

/  .OR.  (  I  V A L (  1  )  . EQ  .  8  )  )  GOTO  91 

IF(  (  1 V  A  L (  5)  .NE.O)  .AND.  (  I  V A L ( 6 )  . N E  .  0  )  ) GOTO  100 

I F (  (  I VAL ( 2 )  . LE .  3 )  .OR.  (  I  V A L ( 2 )  . G E . 9  )  )  GOTO  92 
I F (  (  I VAL( 2 )  . EQ . 4  ) 

/  .OR.  <  I VAL ( 2  )  .  EQ .  6  ) 

/  .OR.  (  I VAL ( 2  )  . EQ . 8 )  )  GOTO  92 

I F (  (  I VAL ( 5 )  . NE . 0 )  .AND.  (  I  V A L (  7 )  . N E . 0  )  ) G OTO  100 


IF( ( IVAL(3) . LE. 3)  .OR.  (  I  V A L ( 3  )  .  GE  .  9  )  )  GOTO  93 
I F (  (  I  V  A L ( 3  )  . EQ . 4 ) 

. OR .  (  I  V  A  L ( 3 )  . EQ  .  6  ) 

•  OR.  ( IVAL( 3)  . EQ.8)  )  GOTO  93 
I F ( (  I VAL( 6 )  . NE . 0  )  .AND.  (  I  V A L ( 8 )  .  N E  .  0  )  )  G OTO  100 

I F (  (  I VAL( 4 )  . LE . 3 )  .OR.  (  I  V A L ( 4 )  .  G E . 9  )  )  GOTO  94 
I F( < I VAL( 4 ) . EQ . 4 ) 

.OR.  ( IVAL( 4 ) .EQ.6) 

•  OR.  (  I  V AL ( 4 )  . EQ . 8 ) )  GOTO  94 
I F (  ( I VAL( 7 )  . NE . 0 )  .AND.  (  I  V  A L ( 8 )  . N E . 0 )  ) G OTO  100 

GOTO  94 

CALL  PPNT ( I BUF  ,  I  I  , J J  ,  1  5  ) 


RETURN 


PROGRAM:  TDT5.FR 


LANGUAGE:  FORTRAN  V 


WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE 

PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  THIN.FR  TO  TEST  FOR 
D  3  DELETION. 

********************************************************* 


SUBROUTINE  TDT5(II,JJ) 

COMMON  I BUF ( 300) 

COMMON  IBUF2O00) 

COMMON  I  B  U  F  3  ( 300) 

COMMON  IARY( 256) 

COMMON  I  NAME (20) 

COMMON  I F  LAG 
COMMON  1ST 
COMMON  LST 
INTEGER  I  V  A  L (  1 2 ) 

DO  5  1-1,12 

I  V  A  L  (  I  )  -  0 
CONTINUE 

CALL  GPNT (  IBUF3,II-1  ,JJ-1  ,  I VAL(  1  )  ) 

CALL  GPNT(  IBUF3  ,  II-l  ,  JJ  +  1  ,  I  VAL( 2  )  ) 

CALL  GPNT(  I BUF3 ,11  +  1  ,  J  J  -  1  ,IVAL(3)> 

CALL  GPNT(  IBUF3  ,  I  1+1  ,  J  J  +  1  ,  I VAL( 4 )  ) 

CALL  GPNT(  I  BUF3  ,  I  I- 1  , JJ ,  I VAL( 5  )  ) 

CALL  GPNT(  I  BUF3  ,  I  I  , J J+  1  ,  I VAL ( 6 )  ) 

CALL  GPNT(  IBUF3  ,  I  I  .  J J- 1  ,  I  VAL( 7  )  ) 

CALL  GPNT  (  I  BUF3  ,  I  I-r  1  ,  J  J  ,  I  VAL(  8  )  ) 

IVAL(9)*IVAL(3)+IVAL(8)+IVAL(4) 

I  V  A  L (  1 0 ) -  I VAL (  1  )  +  I  V  A  L  ( 7 )  +  IVAL( 3  ) 

I  V  A  L (  1  1  )»IVAL(2)  +  IVAL(6)  +  IVAL(4) 

I  V  A  L (  1 2 ) -  I VAL<  1  )  +  IVAL(5)  +  IVAL(2) 

DO  10  1-1,4 

I F ( I VAL( I ) . NE . 8 )GOTO  10 

I  F  (  (  I  VAL  (  I +4  )  .  NE  .  0  )  .AND.  (  I  V  A  L  (  I +8  )  .  E  0 . 0  )  .)  GOTO  100 
CONTINUE 

GOTO  91 

CALL  P PNT (  I BUF  ,  I  I  , J J  ,  l  5  ) 


RETURN 


£***************************************************************£ 

C  C 

C  PROGRAM:  TDT6.FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  THIN. F  R  TO  TEST  FOR  C 
C  D4  DELETION.  C 

C  C 

£***************************************************************£ 
C  C 

SUBROUTINE  TDT6(II,JJ) 

C  C 


COMMON  IBUF( 300) 

COMMON  IBUF2( 300) 

COMMON  IBUF3( 300) 

COMMON  IARY(256) 

COMMON  I  NAME (20) 

COMMON  I F  LAG 
COMMON  1ST 
COMMON  LST 
INTEGER  I  V A L (  1  2  ) 

C 

DO  5  1-1,12 

I  V A L (  I  ) -0 
5  CONTINUE 

CALL  GPNT(  IBUF3  ,  II-l  ,JJ-1  ,IVAL(  1  )) 

CALL  GPNT (  IBUF3  ,  I  I  —  1  , JJ+1  ,  I VAL( 2 )  ) 

CALL  GPNT(IBUF3,II+1  ,  J  J  -  1  ,  I  VAL(  3  )  ) 

CALL  GPNT(IBUF3,  11+ l  , JJ+1  ,IVAL(4)> 

CALL  GPNT(  I BUF3  ,  I  I  -  1  ,  J J  ,  I VAL(  5 )  ) 

CALL  GPNT (  IBUF3,  II  ,  J  J -  l  ,  I  VAL( 6 )  ) 

CALL  GPNT ( I BUF3  ,  I  I  , J J+  1  ,  I VAL( 7  )  ) 

CALL  GPNT(  I BUF3  ,  I  1+  1  ,  J J  ,  I  VAL( 8 ) ) 

IVAL( 9)=IVAL( 7 )+IVAL( 8 )+IVAL( 4) 

IVAL( 10)=IVAL(6)+IVAL(8)+IVAL( 3) 

I  V  A  L (  1  1  )  =  IVAL( 5 )  +  IVAL( 7 )  +  IVAL( 2 ) 

I  V  A  L (  1 2 ) *  I VAL (  1  )  +  IVAL(5)  +  IVAL(6) 

I F (  I VA  L (  1  )  . LE .  5 )  GOTO  91 
IF(  I  V  A  L ( 9 )  .NE.O)  GOTO  91 

I F (  (  I VAL ( 5 )  . NE  .  0  )  .AND.  (  I  V A L ( 6  )  . N E  .  0  )  )  G OT 0  100 

91  I F(  I VAL( 2  )  . LE .  5 )  GOTO  92 
IF( IVAL( 10) .NE.O)  GOTO  92 

IF(  (  IVAL( 5)  .NE.O)  .AND.  (  I  V A L ( 7 )  .  N E  .  0  )  )  G OT 0  100 

92  IF(  IVAL( 3)  . LE.  5  )  GOTO  93 
IF(  IVAL(  1  1  )  .NE.O)  GOTO  93 

I F (  (  I  VAL( 6  )  . NE  .  0  )  .AND.  (  I  V A L (  8 )  . N E .  0  )  ) G OT 0  100 


1  2  1 


I F ( I VAL( 4 ) . LE . 5 )  GOTO  94 
I F ( I VAL (  1 2 )  . NE . 0  )  GOTO  94 

IF( ( I VAL( 7 )  . NE . 0 )  .AND.  (  I  V A L ( 8  )  . N E  .  0  )  ) G OTO  100 
GOTO  94 

CALL  PPNT (  I BUF  ,  I  I  , J J  ,  1  5  ) 

RETURN 

END 


£***********★******************★***■*****************************£ 

C  C 

c  c 

C  PROGRAM:  TGEN.FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  PROGRAM  WILL  COPY  A  CHARACTER  FROM  TSET2.VD  C 

C  (SPECIFIED  BY  THE  USER)  ONTO  A  BLANK  GACKGROUND.  C 

C  THE  USER  MAY  REPEAT  THE  PROCESS  TO  GENERATE  A  C 

C  SEQUENCE  OF  ISOLATED  LETTERS  OR  CONCATENATED  C 

C  LETTERS.  C 

C  C 

C  C 

£***************************************************************Q 

C  C 

c  c 


PARAMETER  I  NAME* "TSET2 . VD" 

PARAMETER  WFILE»"WT" 

IMPLICIT  INTEGER(C,D,E) 

INTEGER  I  COL (  1 6 )  ,  I  LEN (  1  6  ) 

INTEGER  IARY(256)  ,  I BUF (  300  )  ,  IBUF2(  300) 
INTEGER  I BR Y (  2  5  6  ) 

C 


C . C 

C  ;  C 

TYPE"PMOVE  WORKING"  ;  I/O  C 

C  ;  C 

C . ; . C 

C  ;  C 

CALL  P ICFMT( I BUF , INAME , IFLG)  ;  C 

CALL  GFMT(  I BUF , NROWS , NCOLS  , NB I TS  ,  I  MODE  )  ;  C 

;  C 

CALL  MAKB(IBUF)  ; P I C B U F  INPUT  C 

CALL  P I C I N (  I BUF ,  I  NAME ,  I FLG )  ;  BLOCK  C 

;  C 

CALL  P ICFMT (  I  BUF2  ,  WF I LE  ,  I  FLG )  ;  C 

CALL  GFMT(IBUF2, NROWS, NCOLS, NBITS, IMODE) 


CALL  M A  K  B (  I  BUF2  )  ;  C 

CALL  PICIN(  IBUF2  ,  WFILE  ,  IFLG)  ;  C 

C  ;  C 

c . ; . c 

C  ;  C 

1  TYPE "ENTER  CHARACTER  TO  MOVE:"  ;  C 

TYPE" 1 ->A ,  2  - >  B  ,  3 - >  C  ,  4  - >  D "  ;  ENTER  LETTER  C 

TYPE "UP  TO  16->P"  ;  C 

ACCEPT" ENTER  NUMBER:  ",INUM  ;  C 

C . ; . C 

C  ;  C 

;  C 

DATA  (  I  CO  L (  I  )  ,  I  *  l  ,  l  6  )  ;  C 

/  /5,«0  ,  7  1  ,  1  0  1  ,  1  3  5  ,  1  6  3,  1  8  9,2  2  1  ,  ;  LETTER  POSITION  C 

/  h ,  2  5  ,  5  3 , 38  ,  1  20  ,  1  60 ,  1  96  ,  2  29/  ;  C 

;  AND  C 


1  2  3 


n  n 


DATA  (  I LEN (  I  )  ,  I  *  1  ,  1  6  )  ;  C 

/  /  2  7 , 22  ,  22  ,  2  5  ,  20  ,  1  9  ,  26 , 2  7  ,  ;  LENGTH  C 

/  14, 16,26,22,30,26,25,21/  ;  C 

C . ; . C 

C  ;  C 

I  TO  P  *  4  2  ;  C 

IF(  INUM.GE.9)IT0P»1  1  3  ;  SET  STARTING  C 

I BOT “63  ;  ROWS  C 

I  F  ( INUM.GE. 9) I BOT= 1 42  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

TYPE  ;  ENTER  C 

TYPE "ENTER  POSITION  TO  PLACE  CHARACTER:"  ;  NEW  POSITON  C 

ACCEPT" ENTER  TOP  ROW  :  " , JTOP  ;  OF  WORD  C 

ACCEPT" ENTER  LEFT  COL:  ",JCOL  ;  C 

C . ; . C 

C  ;  C 

;  C 

I LN= I LEN ( INUM)  ;  MOVE  LETTER  C 

IC*ICOL( INUM)  ;  TO  C 

J-JTOP  ;  C 

DO  10  I s ITOP , I BOT  ;  NEW  C 

CALL  CROW(  IBUF  ,  I  ,  IC ,  I LN ,  IARY)  ;  POSITION  C 

CALL  GROW(  IBUF2  ,J , JCOL ,  ILN  ,  I BRY )  ;  C 

;  C 

;  C 

DO  15  K* 1 , I LEN  ;  C 

I F (  (  IARY(K)  .  EQ.  0)  ;  C 

/  .  OR .  ;  C 

/  (  IBRY(K)  . EQ . 0 ) )IBRY(K)*0;  C 

15  CONTINUE  ;  C 

;  c 

CALL  PROW(  I  BUF 2  ,  J  ,  JCOL  ,  I  LN  ,  I  BRY  )  ;  C 

J-J+l  ;  C 

1  C 

10  CONTINUE  ;  C 

C  ;  C 

C . ; . C 

C  ;  C 

TYPE  ;  C 

ACCEPT "CONTINUE  (Y»l):  ",ICNT  ;  REPEAT  FOR  MORE  C 

IF(  ICNT.EQ.  1  )GOTO  1  ;  LETTERS  C 

C  ;  C 

C . ; . C 

C  ;  C 

CALL  PICOUT( I  B  U  F  2 , WFI LE , IFLG)  ;  PICBUF  OUTPUT  C 

CALL  RELB(IBUF)  ;  BLOCK  C 

CALL  RELB(  IBUF2  )  ;  C 

C  ;  C 

C . ; . C 

C  C 

STOP 
END 


1  2  4 


£***************************************************************£ 

c  c 

c  c 

C  PROGRAM:  TT2.FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN :  1 5  DEC  84  SYSTEM:  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  THIN. FR  TO  TEST  FOR  C 
C  P I  DELETION  .  C 

C  C 

C  C 

^  *  *********■**•******************★*****★****★*★.  *****★******★*****£ 

C  C 

SUBROUTINE  TT2(II,JJ) 

C  C 

C  C 


COMMON  I  BUF (  300  ) 

COMMON  I BUF2 ( 300 ) 

COMMON  IBUF3(300) 

COMMON  IaRY(256) 

COMMON  I  NAME (20) 

COMMON  I F  LAG 
COMMON  1ST 
COMMON  LST 
INTEGER  I  V  A  L ( 9  ) 

C 

C 

DO  5  1-1,9 

I  V  A  L (  I  )  -0 

C  •  5  CONTINUE 

CALL  GPNT(  IBUF3  ,  I  I  —  I  , J  J -  1  ,  IVAL(  1  )) 
CALL  GPNT(IBUF3,II-1,JJ,IVAL(2)) 
CALL  GPNT(  IBUF3  ,  I  I -  1  , JJ+1  ,  I VAL (  3 )  ) 

CALL  G  P  N  T  (  IBUF3,  I  I  ,  JJ  -  1  ,  I  VAL(  *.  )  ) 
CALL  GPNT (  IBUF3,  II  , J  J+ 1  ,  I  V  A  L (  5 )  ) 

CALL  G  P  N  T (  I  BUF  3  ,  I  I  +  1  ,JJ-1  ,  I  VAL (  6  )  ) 

CALL  G PNT (  I B U F 3  ,  I  I  +  1  , J J  ,  I  V A L  (  7  )  ) 

CALL  G  P  NT (  IBUF3,  I  1  +  1  ,  J  J  +  1  ,  I  V AL ( 8  )  ) 


DO  10  1=1,8 

I  F(  I  VAL(  I  )  . EQ .  2  )GOTO  100 
I F (  I  V  A  L (  I  )  .EQ  5  ) G 0 T  0  10  0 

I F ( I VAL ( I ) . EQ . 6 ) GOTO  100 
IF(IVAL(I).EQ.7)C0T0  100 
10  CONTINUE 

GOTO  94 

100  CALL  PPNT(  IBUF,  I  I  , J J  ,  13) 

9  h  RETURN 

END 


1  2  5 


(-*************************************************************«*£ 

c  c 

C  PROGRAM:  TT3.FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN:  15  DEC  8-.  SYSTEM  :  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  THIN.FR  TO  TEST  FOR  C 
C  P 1  DELETION  .  C 

C  C 

C  C 

£***************************************************************£ 

c  c 

SUBROUTINE  TT3(II,JJ) 

C  C 


COMMON  I  BUF (  300  ) 

COMMON  IBUF2(300) 

COMMON  I  BUF3  (  300  ) 

COMMON  I  ARY (  2  5  6  ) 

COMMON  I  NAME ( 20 ) 

COMMON  I  FLAG 
COMMON  1ST 
COMMON  LST 
INTEGER  I  V  A  L (  9  ) 

C 

C 

DO  5  1*1,9 

I  V  A  L  (  I  )  *  0 
5  CONTINUE 

CALL  GPNT (  IBUF3,  I  I  —  1  , J J  ,  I VAL (  1  )  ) 

CALL  G  P  N  T (  I8UF3,  I  I  ,  J  J  -  1  ,  I  V  A  L (  2  )  ) 

CALL  GPNT (  IBUF3  ,  I  I  , JJ+l  ,  I  V  A  L (  3  )  ) 

CALL  GPNT(  IBUF3,  II+l  ,JJ,  I  V  A  L ( 4  )  ) 

DO  10  I « 1 , h 

IF( I VAL ( I ) , EQ . m ) GOTO  90 
10  CONTINUE 

GOTO  9x 


9  0 

I  F  (  (  I  V  A  L  ( 

GOTO  92 

1  )  .  EQ .  a  ) 

.  OR  . 

(  I  VAL (  a )  .  EQ .  4 )  ) GOTO 

9  1 

9  1 

IF(  (  I  V  A  L ( 

2 ) . NE . 0 ) 

.OR. 

(  I  VAL (  3)  . NE .  0 )  ) GOTO 

1  00 

9  2 

I F (  (  I  VAL (  2 )  ,  EO  .  4  ) 

GOTO  9a 

.  OR  . 

(  I  V  A  L  (  3  )  .  E  Q  .  4  )  )  G  0  T  0 

9  3 

9  3 

I  F  (  (  I  V  A  L  ( 

1  )  .  N  E  .  0  ) 

.  OR  . 

(  I  V  A  L ( 4 )  .NE.O)  ) GOTO 

1  00 

GOTO  9a 


100  CALL  P  P NT (  I  B  U  F ,  I  I  ,  J  J  ,  1  5  ) 

9a  RETURN 

END 
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I  S 


^  in 
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;  !  :  :  -  r  i  ,  t  t  h  e  -  n  g  1  i  s  h  alphabet.  The  letters  :n  a  v  be 

isolated  letters;  or  concatenated,  groups  of  letters 
s  1  g  >  r  i  _  h  si  essentially  consists  of  a  two  pass  i  o  r  r  e  1  a  t 

w  h  1  h  ;  o  .»  s  :  e  in  d  e  n  e  r  g  y  of  the  unknown  chariot  •  r  s  at 
j  i  a  1  to  a  known  template. 
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ion  in 

<i  set 

i  a 
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Lugene  Shura  was  born  in  New  York.  City,  New  York,  on  20 
April  196L.  He  graduated  from  Manhattan  College,  New  York, 
in  1983.  Upon  graduation,  he  was  commissioned  as  a  Second 
Lieutenant  th  ough  the  AFROTC  program  and  was  assigned  to 
the  School  of  Engineering,  Air  Force  Institute  of 
Technology.  Upon  graduation  from  AFIT,  he  was  assigned  to 
Edwards  Air  Force  Base. 

Permanent  Address:  50  Bayard  St. 

Apt.  2  H 

New  York,  New  York  10013 


0* 


TYPE"VTHIN 

WORKING"  ; 

I/O 

C 

C 

DO  1  ICNT-1 

,  MAX 

REPEAT 

c 

PROCESS  MAX 

c 

TIMES 

c 

CALL  MTRXA 

FORM 

c 

MATRIX  A 

c 

DO  5  I I - 1ST 

,  LST 

c 

c 

DO  10  JJ 

=*2,255 

c 

CALL 

G  PNT ( I BUF3 , I I  , J J ,  I V ) 

c 

GOTO 

(81,82,83,84,85,86,87,88,  10)IV 

c 

c 

GOTO 

1  0 

c 

c 

8  I 

CALL 

P  P  NT (  I BUF ,  I I  , J J  ,  1  5  ) 

FILL  IN 

c 

GOTO 

1  0 

HOLES 

c 

c 

82 

CALL 

TT2 ( I I , J J ) 

TEST  P  2 

c 

GOTO 

1  0 

c 

c 

83 

CALL 

TT  3 (  I I  , J  J  ) 

TEST  P  3 

c 

CALL 

TDT3 ( I I , J J ) 

TEST  D  3 

c 

GOTO 

l  0 

c 

c 

84 

CALL 

TT  4 (  I I  , J  J ) 

TEST  P  4 

c 

CALL 

TDT  4 ( I I  , J J ) 

TEST  D  4 

c 

GOTO 

10 

c 

c 

8  5 

CALL 

TT5( II , JJ) 

TEST  P  5 

c 

CALL 

TDT  5 ( I I , J J ) 

TEST  D  5 

c 

GOTO 

10 

c 

c 

86 

CALL 

TT 6 (  I I  , JJ  ) 

TEST  P  6 

c 

CALL 

TDT 6 (  I I  , J J  ) 

TEST  D  6 

c 

GOTO 

10 

c 

c 

87 

CALL 

TT  7 ( I I  , JJ) 

TEST  P  7 

c 

GOTO 

10 

c 

1 

c 

88 

CALL 

TT8 ( I I , J J ) 

TEST  P  8 

c 

1  0 

CONTINUE 

c 

c 

5 

CONTINUE 

c 

I 

CONTINUE 

c 

c 

9999 

CALL  PICOUT(  IBUF,  INAME,  IFLG) 

* 

;  PICBUF  OUTPL 

CALL  RELB(IBUF) 

9 

BLOCK 

CALL  RELB(IBUF2) 

CALL  RELB( I B UF  3  ) 

C . 

c 

ST0P"<7>THIN  FINISHED" 
END 


1 

I 

J 
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************«**************************************************(-- 

PROGRAM:  VTHIN.FR  LANGUAGE:  FORTRAN  V  C 


LANGUAGE:  FORTRAN  V 


WRITTEN:  15  DEC  84 


SYSTEM:  ECLIPSE 


C  PURPOSE:  THIS  PROGRAM  THINS  A  PATTERN  TO  PRODUCE  A  SKELETON. C 
C  C 

C  NOTE : 1 .  THE  INPUT  PATTERN  MUST  BE  BINARY  (0/15  FOR  C 

C  THE  OCTEK ) .  C 

C  2.  THE  INPUT  PATTERN  MUST  FIRST  BE  PRE-PROCESSED  C 

C  USING  SPOTS.FR  AND  MEDIAN. FR.  C 

C  3.  SOURCE  OF  ALGORITHM:  C 

C  LINE  THINNING  ALGORITHM  C 

C  G.  FEIGIN  AND  N.  BEN-YOSEF  C 

C  PROCEEDINGS  OF  SPIE,  VOL  397  C 

£***************************************************************£ 
C  C 

COMMON  I  BUF(  300  ) 

COMMON  I  BUF2  (  300  ) 

COMMON  I  BUF  3 (  300  ) 

COMMON  I AR  Y (  2  5  6  ) 

COMMON  I  NAME (  20 ) 

COMMON  I F  LAG 
COMMON  1ST 
COMMON  LST 
REAL  W  F I L  E (  5  ) 


INTEGER  MAX 


ACCEPT" ENTER  INPUT  FILENAME:  " 
R  E  A  D (  11,199)  I N  AME (  1  ) 

FORMAT ( S  *0 ) 


ENTER 
F I LENAME 


ACCEPT" ENTER  STARTING  ROW:  ".1ST 
ACCEPT "ENTER  END  ROW:  ",LST 


ENTER  ROW 
POSITIONS 


;  ENTER 
;  THICKNESS 

ACCEPT "ENTER  NUMBER  OF  ITERATIONS:  ".MAX  ;  OF  CHARACTER 


WFILE»"BK" 

CALL  P  ICFMT (  I BUF  ,  INAME ,  I FLG ) 

CALL  GFMT( IBUF ,NROWS , NCOLS .NBITS , IMODE) 
CALL  MAKB ( IBUF) 

CALL  P  I  C I N (  I  BUF  ,  I  N AME ,  I FLG ) 

CALL  P  ICFMT(  I BUF2 , WF I LE ,  I FLG ) 

CALL  G  FMT (  IB UF2.NR0WS, NCOLS  .NBITS,  I  MODE) 
CALL  MAK  B (  IBUF2 ) 

CALL  P  IC I N (  I BUF2 , WF I LE ,  I FLG ) 

CALL  PICFMT( IBUF3.WFILE, IFLG) 

CALL  GFMT(IBUF3,NR0WS, NCOLS, NBITS,  I  MODE) 
CALL  MAK  B (  I BUF  3 ) 

CALL  PICIN( IBUF3 ,WFILE, IFLG) 


PICBUF  INPUT 
BLOCK 


c***************************************************************c 

c  c 

C  ORIGINAL  PROGRAM:  VTOC.FR  BY  R.  SIMMONS  C 

C  ORIGINAL  PROGRAM:  UNPACK. FR  BY  R.  SIMMONS  C 


C  PROGRAM:  VTC.FR 


LANGUAGE:  FORTRAN  V 


C  WRITTEN:  15  DEC  84 


SYSTEM:  -CLIPSE 


C  PURPOSE:  THIS  PROGRAM  CONVERTS  A  256X256  PACKED  VIDEO  FILE  C 
C  INTO  A  256X256  UNPACKED  COMPLEX  VIDEO  FILE.  C 

C  C 

c***************************************************************c 

c  c 

INTEGER  I FLNM( 7 ) , OFLNM( 7 ) 

INTEGER  TEMP ( 5 1 2 ) , I N ( 2 048 ) 

COMPLEX  OUT(2048) 


ACCEPT  "ENTER  INPUT  FILENAME:  " 
READ  (11,199)  I  F  LNM (  1  ) 

ACCEPT  "ENTER  OUTPUT  FILENAME: 
READ  (11, 1 99 ) OFLNM( 1 ) 

FORMAT  ( S  4  0 ) 


ENTER 

INPUT  AND  OUTPUT 
F I LENAMES 


OPEN  1 , IFLNM 

CALL  DFILW  (OFLNM.IER) 

I F ( IER.NE. 1 ) GOTO  90 
CALL  CFI LW( OFLNM , 2 , KER > 

CALL  CHECK(KER) 

OPEN  2  ,  OFLNM , ATT* "OR"  , ERR- 1 00 


OPEN 

INPUT  AND  OUTPUT 
FLIES 


DO  5  1-0,31 

CALL  RDBLK( 1 ,2*1, TEMP, 2 ,IER) 

CALL  CHECK ( I ER ) 

CALL  UNPACK( 5 1 2 ,TEMP , IN) 

DO  3  J- 1 , 2048 

OUT(J)-CMPLX(FLOAT( IN(J) ) ,0.0) 
CONTINUE 

CALL  WRBLK(2  ,  3  2*1  .OUT,  32  ,  IER) 

CALL  CHECK (IER) 

CONTINUE 


CONVERT  INTEGER 
INTO  COMPLEX 
(A  +  j  0.0) 


CALL  RESET 

STOP  " <  7  >  FINISHED" 

END 
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PROGRAM:  TT8 . FR 


LANGUAGE:  FORTRAN  V 


WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE 

PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  THIN.FR  TO  TEST  FOR 
P  8  DELETION. 

***************************************************** 


SUBROUTINE  TT8(II,JJ) 


C 


COMMON  IBUF(300) 

COMMON  IBUF2( 300) 

COMMON  IBUF3(300) 

COMMON  I  ARY (  256  ) 

COMMON  I  NAME (20) 

COMMON  I  FLAG 
COMMON  1ST 
COMMON  LST 
INTEGER  I  V A L ( 9 ) 

C 


DO  5  I  ■*  l  ,  9 

I  VAL  (  I  )  =* 0 
CONTINUE 


CALL  GPNT( IBUF3  ,  I  I  —  1  ,  J  J  ,  IVAL(  1  )  ) 
CALL  GPNT(  IBUF3  ,  II  ,  J  J -  1  ,  I  V  A  L ( 2 )  ) 
CALL  GPNT( I BUF3  ,  I  I  , JJ+1  ,  I VAL( 3 )  ) 
CALL  GPNT(  I BUF3  ,  I  1  +  1  ,  J J  ,  I  VAL( 4  )  ) 

CALL  GPNT(IBUF3,II-1  ,JJ-1  ,  IVAL( 5  >  ) 
CALL  GPNT(  IBUF3,  I  I  —  1  ,  J  J+l  ,  I VAL( 6 )  ) 
CALL  GPNT( IBUF3 ,11  +  1  ,JJ-1  ,  I VAL( 7 )  ) 
CALL  GPNT(  IBUF3  ,  I  I  +  1  ,  J  J  +  1  ,  I VAL( 8 )  ) 


I F ( I VAL( 1 ) . NE . 7 )  GOTO  91 
IF( IVAL( 2 ) . NE . 7 )  GOTO  91 

I F(  ( I VAL( 6 )  . EQ . 9 )  .OR.  ( I VAL ( 3 )  . EQ . 9 )  .OR.  ( I  V AL ( 8 )  . EQ . 9 ) 
.OR.  (  IVAL( 4)  . EQ . 9 )  .OR.  (  I  V A L ( 7  )  .  E 0 . 9  )  ) G OTO  100 


IF ( I VAL( 1 ) . NE . 7 )  GOTO  92 
IF(  I VAL( 3 )  . NE . 7  )  GOTO  92 

I F (  (  I VAL ( 5 )  . EQ . 9 )  . OR . ( I VAL ( 2  )  . EQ . 9  )  .OR.  ( I  V A L ( 7  )  . E Q . 9  ) 
.OR.  (  I VAL( 4 )  . EQ . 9 )  . OR .  ( I VAL( 8 )  . EQ  .  9  )  )GOTO  100 


I F ( I VAL( 4 ) . NE . 7 )  GOTO  93 
I F ( I VAL( 2 ) . NE . 7 )  GOTO  93 

IF(  (  I VAL( 5 )  . EQ . 9 )  . OR .  ( I VAL(  1  )  . EQ . 9  )  .OR.  ( I  V A L ( 6 )  . E Q . 9 ) 
.OR.  (  I  V  A  L ( 3 )  .EQ.9)  .OR.  (  I  V  A  L ( 8 )  . EQ . 9  )  )GOTO  100 


I  F  (  I VAL ( 4 )  . NE . 7 )  GOTO  94 
I F ( I VAL( 3 ) . NE . 7 )  GOTO  94 

I F (  (  I VAL( 7 )  . EQ . 9 )  . OR .  ( I VAL( 2  )  . EQ . 9  )  .OR.  (  I  V A L ( 5  )  . EQ . 9  ) 
.OR.  (  I  V  A  L (  1  )  .EQ.9)  .OR.  (  I  V  A  L ( 6 )  .EQ.9))GOTO  100 


X F ( I VAl( 4 ) . NE . 9 )  GOTO  94 
IF(IVAL( 1 ) . NE . 0 )  GOTO  94 

IF(((IVAL(6).EQ.O).OR.(IVAL(3).NE.O))  .AND. 
<(IVAL(5).EQ.O).OR.(IVAL(2).NE.O))  )GOTO  1 0  0 

I F ( I VAL( 2 ) . NE . 8 )  GOTO  95 
IF( I VAL( 7 )  . NE . 7  )  GOTO  95 

I F (  (  I  V  A  L ( 5 )  . EQ . 9 )  .OR.  (  IVAL(  1  )  .EQ.9)  )GOTO  100 

IF( IVAL( 3) . NE. 8)  GOTO  96 
I F ( I VAL( 8 ) . NE . 7 )  GOTO  96 

IF(  ( IVAL(  1  )  . EQ. 9)  .OR .  ( IVAL(6)  . EQ. 9)  )GOTO  100 

IF( IVAL(2) .NE.8)  GOTO  97 
IF( IVAL( 5 ) . NE . 7 )  GOTO  97 

I F ( (  I  V  A  L ( 7 )  .EQ.9)  .OR.  (  IVAL(4)  .EQ.9))G0T0  100 


IF(IVAL(3) .NE.8)  GOTO  98 
I F ( I VAL( 6 ) . NE . 7 )  GOTO  98 

IF((IVAL(4). EQ.9). OR.  (IVAL(8). EQ.9)  )GOTO  1 0  0 


GOTO  98 


CALL  PPNT( I BUF  .  I  I  ,  J J  ,  1  5  ) 

RETURN 

END 


PROGRAM:  TT7.FR 


LANGUAGE:  FORTRAN  V 


WRITTEN:  15  DEC  84 


SYSTEM:  ECLIPSE 


PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  THIN.FR  TO  TEST  FOR 
P  8  DELETION. 


SUBROUTINE  TT7(II,JJ) 


COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

COMMON 

INTEGER 


I  BUF (  300) 
I8UF2(300) 
IBUF3(300) 
IARY(  2  5  6  ) 

I  NAME ( 20 ) 

I  F  LAG 
I  ST 
LST 

:  I  V  A  L  (  9  ) 


DO  5  1-1,9 

I  V  A  L  (  I  )  -  0 
CONTINUE 

CALL  GPNT( I BUF3 ,  I  I -  1  , JJ  ,  I VAL(  1  ) ) 

CALL  GPNT(  IBUF3 ,  II  ,  J  J -  1  ,  I  V  AL ( 2 )  ) 

CALL  GPNT(  I  BUF3 ,  I  I  ,  J J+l  ,  I  VAL( 3 )  ) 

CALL  GPNT(  I  BUF3 ,  I  1  +  1  , JJ ,  I VAL( 4 )  ) 

CALL  GPNT(  IBUF3,  II-l  , JJ-1  ,  I VAL( 5 )  ) 

CALL  GPNT (  IBUF3 ,  II-l  ,  J  J+  1  ,  IVAL(6)) 

CALL  GPNT(  IBUF3  ,  II  +  l  ,  JJ-1  ,  I VAL( 7 )  ) 

CALL  GPNT(  IBUF3 ,  II+l  , JJ+1  ,  I VAL( 8 )  ) 

I  F (  I  VAL(  l  )  . NE . 9 )  GOTO  91 
IF(  I  VAL( 4 )  .NE.O)  GOTO  91 

I  F ( (  ( I  V  A  L ( 7 )  .EQ.O)  .OR.  ( I  V  A  L ( 2 )  .NE.O)  )  .AND. 
(  (  IVAL(8)  .  EQ.O)  .OR.  (  I  V  A  L ( 3 )  .NE.O)  )  )GOTO  100 


IF (  I  VAL( 2 )  . NE . 9 )  GOTO  92 
I  F  (  I VAL( 3 )  . NE . 0 )  GOTO  92 

IF(((IVAL(8). EQ.O). OR. ( I VAL( 4 ) . NE.O))  .AND. 
((IVAL(6).EQ.O).OR.(IVAL(L).NE.O)))GOTO  100 


I F  (  I VAL  (  3 )  . NE  .  9 )  GOTO  93 
IF( I  V  A  L ( 2) .NE.O)  GOTO  93 

IF(((IVAL(5). EQ.O). OR. (IVAL(l). NE.O))  .AND. 
(  (  I  V  A  L ( 7 )  .EQ.O)  .OR.  (  I  V  A  L ( 4 )  .NE.O)  )  )GOTO  100 


PROGRAM:  TT6 . PR 


LANGUAGE:  FORTRAN  V 


WRITTEN:  15  DEC  8  4  SYSTEM:  ECLIPSE 

PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  THIN.FR  TO  TEST  FOR 
P  7  AND  NEIGHBOR  VALUE  OF  9  DELETION. 

ft********************************************************** 


SUBROUTINE  TT6(II,JJ) 

COMMON  IBUF(300) 

COMMON  I BUF2 (  300  ) 

COMMON  IBUF3C300) 

COMMON  IARY(256) 

COMMON  INAME( 20  ) 

COMMON  I  FLAG 
COMMON  1ST 
COMMON  LST 
INTEGER  I  V A  L ( 9  ) 

INTEGER  T ( 3  ) 

DO  5  1*1,9 

I  V  A  L (  I  )-0 
CONTINUE 

CALL  GPNT(  IBUF3 ,  I  I  —  1  , J J , IVAL(  1  )  ) 

CALL  GPNT(  IBUF3  ,  II  ,  J  J -  1  ,  I  V  A  L ( 2  )  ) 

CALL  GPNT( IBUF3  ,  I  I  , JJ+1  ,  IVAL( 3  )  ) 

CALL  GPNT(  I BUF3  ,  I  1  +  1  ,  J J  ,  I  VAL( 4  )  ) 

DO  7  1-1,4 

IF (  I  VAL(  I )  . EQ . 9 )GOTO  100 
CONTINUE 

CALL  GROW(  I BUF3 ,  I  1  +  1  , J J- 1  , 3 , T) 

IVAL(5)-T( 1 )+T(2)+T(3) 

CALL  GROW( I BUF3  ,  I  I -  1  ,  J J-  1  , 3 , T ) 

IVAL(8)-T(1)+T(2)+T(3) 

CALL  GPNT(  IBUF3  ,  I  I  -  1  , JJ+l  ,T(  1  ) ) 

CALL  GPNT( I BUF3  ,  I  1  +  1  , J J+l  , T< 2 )  ) 

IVAL(6)-T( 1 ) +T (2)+IVAL(3) 

CALL  GPNT( IBUF3  ,  I  I  —  1  ,  J  J  -  1  ,T(  1  )  ) 

CALL  GPNT(  IBUF3,  I  I  +  1  ,  J  J  -  1  ,T( 2)  ) 
IVAL(7)-T(1)+T(2)+IVAL(2) 

DO  10  1-1,4 

IF( I VAL( I  )  . LE .  5 )  GOTO  10 

IF(  (  I VAL(  I )  . EQ .  6 )  .OR.  (  I  V A L (  I )  . EQ . 8  )  )  GOTO  10 
I F(  I VAL( 1  +  4 )  . EQ  .  0  )  GOTO  100 
CONTINUE 
GOTO  91 

0  CALL  PPNT(  IBUF ,  I  I  ,  JJ  ,  1  5  ) 

RETURN 

END 


0 


3 


IF( IVAL( 3 ) . NE . 5 )GOTO  83 

IF( ( I VAL(  1  ) . EQ. 0 )  .OR.  (  I  V A L ( 5 ) . E Q . 0 ) ) GOTO  83 
IF( I VAL( 7 ) . EQ . 0 )GOTO  100 

I F( I VAL( 4 ) . NE . 5 )GOTO  84 

I F ( ( I VAL( 3 )  . EQ . 0 )  .OR.  ( I  V A L ( 6 )  . EQ . 0 ) ) G OTO  84 
I F ( I VAL( 5 ) . EQ . 0 )GOTO  100 

IVAL(9)-IVAL(4)+IVAL(7)+IVAL(8) 

IVAL(10)-IVAL(3)+IVAL(6)+IVAL(8) 

IVAL(11)»IVAL(2)+IVAL(5)+IVAL(7) 

IVAL( 12)-IVAL( 1 )  +  I  V  A  L  ( 5)+IVAL(6) 


DO  10  1-1,4 

IF((IVAL(I) . LT  .  5  )  .OR.  ( I VAL ( I ) . GT . 7 ) ) GOTO  10 
I F( I VAL( 1+8 ) . EQ . 0 )GOTO  100 
CONTINUE 


GOTO  91 

CALL  PPNT( I BUF  ,  I  I  ,  J J  ,  l  5  ) 


RETURN 


PROGRAM:  TT5.FR 


LANGUAGE:  FORTRAN  V 


WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE 


PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  THIN.FR  TO  TEST  FOR 
P  5  ,  P  6 ,  AND  NEIGHBOR  VALUE  OF  8  DELETION. 


********************************************************* 

SUBROUTINE  TT5(II,JJ) 

COMMON  IBUF(300) 

COMMON  I BUF2 (  300  ) 

COMMON  IBUF3( 300) 

COMMON  IARY(256) 

COMMON  INAME(20) 

COMMON  I F  LAG 
COMMON  1ST 
COMMON  LST 
INTEGER  I  V A L (  1  2  ) 

DO  5  1*1,9 

I  V  A  L  (  I  )  “  0 
CONT I NUE 

CALL  GPNT (  I BUF3 ,  I  I - 1  , J J ,  I VAL(  1  )  ) 

CALL  GPNT(  IBUF3,  II  ,  J  J - l  ,  I  V  A  L ( 2 )  ) 

CALL  GPNT ( IBUF3 ,  I  I , J J+l  ,IVAL(3)) 

CALL  GPNT(  IBUF3,  II+l  ,  J  J  ,  I  V  AL ( 4 )  ) 


DO  7  1-1,4 

IF( IVAL( I ) . EQ . 8)  GOTO  100 
CONTINUE 


CALL  GPNT (  IBUF3,  I  I  —  1  ,JJ-1  ,  I VAL(  5  )  ) 
CALL  GPNT(  IBUF3  ,  I  1-1  ,  J  J  +  1  ,  IVAL( 6)  ) 
CALL  GPNT ( IBUF3 , II+l  ,  J J- 1  , I VAL( 7  )  ) 
CALL  GPNT(  IBUF3  ,  II  +  l  , JJ+1  ,  I VAL( 8  )  ) 


IF( IVAL( 1 ) . NE . 5 )GOTO  81 

I F (  (  I VAL( 2 )  . EQ . 0 )  .OR.  (  I  V A L ( 7 )  . EO  .  0  )  )  GOTO  81 
IF( I VAL( 8 ) . EQ . 0 )GOTO  100 


I F ( I VAL( 2 ) . NE . 5 )GOTO  82 

I F (  (  I VAL( 4 )  . EQ . 0 )  .OR.  (  I  V A L ( 8  )  . EQ  .  0  )  ) GOTO  82 
I F ( I VAL( 6 ) . EQ . 0 )GOTO  100 


£***************************************************************Q 

C  C 

c  c 

C  PROGRAM:  TT4.FR  LANGUAGE:  FORTRAN  V  C 

C  C 

C  WRITTEN:  15  DEC  84  SYSTEM:  ECLIPSE  C 

C  C 

C  PURPOSE:  THIS  IS  A  SUBROUTINE  USED  BY  THIN.FR  TO  TEST  FOR  C 

C  P  2 ,  P  3 ,  P  4  DELETION.  C 

C  C 

C  C 

0***************************************************************C 

c  c 

SUBROUTINE  TT4(II,JJ) 

C  C 


COMMON  IBUF(300) 

COMMON  I8UF2(300) 

COMMON  IBUF3(300) 

COMMON  IARY(256) 

COMMON  I  NAME (20) 

COMMON  I F  LAG 
COMMON  1ST 
COMMON  LST 
INTEGER  I  V A L ( 9 ) 

C 

DO  5  1-1,9 

I VAL (  I  )  -0 
5  CONTINUE 

CALL  GPNT(  IBUF3  ,  I  I  —  1  ,  J  J  ,  I  V  A  L (  1  )  ) 

CALL  GPNT(  IBUF3  ,  I  I  ,  JJ-  1  ,  IVAL( 2 )  ) 

CALL  GPNT(  IBUF3  ,  I  I  , J J+l  ,IVAL(3)) 

CALL  GPNT( IBUF3,  I  I  +  1  ,JJ,  I  V  A  L ( 4 )  ) 

I F (  ( I VAL ( 2 )  . EQ . 0 )  .OR.  (  I  V AL (  3 )  . EQ . 0 ) ) GOTO  91 
IF((IVAL(  l  )  . G  E  .  4  )  .AND.  (  I  V A L (  1  )  . L E .  7 )  ) G OTO  100 

IF(  (  I  V  A  L ( 4 )  . G  E . 4  )  .AND.  ( I VAL ( 4 )  . LE . 7  )  )  GOTO  100 

91  IF( (  IVAL(  1  )  . EQ . 0 )  .OR.  ( I  V A L < 4 )  . EQ . 0  )  ) GOTO  92 

I F (  (  I VAL ( 2 )  . GE . 4  )  .AND.  (  I  V A L ( 2 )  .  L E . 7  )  )  G OTO  100 

IF( (  IVAL( 3)  .GE .  4)  .AND.  (  I  V A L ( 3  )  . L E .  7  )  ) G OTO  100 

92  I F (  (  I  V  A  L ( 2 )  .NE.O)  .OR.  (  I  V A L ( 3  )  .  N E  .  0  )  ) GOT 0  93 

I F (  (  I  V  A  L (  l  )  . EQ  .  7 )  .OR.  (  I  V A L ( 4  )  . EQ  .  7  )  ) GOTO  100 

93  IF(  (  IVAL(  1  )  .NE.O)  .OR.  (  I  V A L ( 4 )  .  N E  .  0  )  ) G OTO  94 

I F (  (  I VAL( 2 )  . EQ .  7 )  .OR.  (  I  V A L ( 3  )  .  EQ .  7  )  )  G OTO  100 
GOTO  94 

100  CALL  PPNT(  IBUF,  II  ,JJ,  1  5) 

94  RETURN 
END 
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